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Abstract 

Motivated by orbifold grand unified theories, we construct a class of three-family Pati-Salam models in 



> 
oo 

On _ 

a Z6 abelian symmetric orbifold with two discrete Wilson lines. These models have marked differences 

ON 

from previously-constructed three-family models in prime-order orbifolds. In the limit where one of the six 
compactificd dimensions (which lies in a Z2 sub-orbifold) is large compared to the string length scale, our 
models reproduce the supersymmetry and gauge symmetry breaking pattern of 5d orbifold grand unified 
1 I theories on an S 1 /Z2 orbicirclc. We find a horizontal 2 + 1 splitting in the chiral matter spectra - 2 families 

of matter are localized on the Z2 orbifold fixed points, and 1 family propagates in the 5d bulk - and identify 



them as the first-two and third families. Remarkably, the first two families enjoy a non-abelian dihedral D4 
family symmetry, due to the geometric setup of the compactificd space. In all our models there are always 



some color triplets, i.e. (6, 1, 1) representations of the Pati-Salam group, survive orbifold projections. They 
could be utilized to spontaneously break the Pati-Salam symmetry to that of the Standard Model. One 
model, with a 5d Eg symmetry, may give rise to interesting low energy phenomenology. We study gauge 
coupling unification, allowed Yukawa couplings and some of their phenomenological consequences. The Eg 
model has a rcnormalizable Yukawa coupling only for the third family. It predicts a gauge- Yukawa unification 
relation at the 5d compactification scale, and is capable of generating reasonable quark/lepton masses and 
mixings. Potential problems are also addressed, they may point to the direction for refining our models. 
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I. INTRODUCTION AND MOTIVATION 



The Standard Model (SM) has been a cornerstone of modern-day particle physics. Although 
during the past three decades it has passed all experimental tests, nevertheless there are many 
open questions remaining to be answered. We have yet to understand — (i) the mechanism of 
electro-weak symmetry breaking, and find the Higgs boson which might be responsible for this 
breaking; (ii) the quantized fermion charges (why the up and down quarks have charges 2/3 and 
—1/3 respectively) and the weak mixing angle (why it is 0.23); (iii) the 3 replicas of quarks and 
leptons, the observed fermion mass hierarchy, and the Cabibbo-Kobayashi-Maskawa (CKM) mixing 
matrix and its leptonic cousin. We ought to understand these problems from more fundamental 
principles, rather than simply take various charges, masses, mixing angles and CP phases as input 
parameters as in the SM. Addressing these questions may eventually lead us to a more fundamental 



theory such as string theory at high energy scales. 

Q 



String theory is a leading candidate for a consistent theory of quantum gravity. It has a rich 
structure and many believe it can easily accommodate the SM as a subset. Moreover there have 
been many attempts in the past to construct supersymmetric generalizations of the SM (which 
will be loosely referred to as the minimal supersymmetric standard model, or the MSSM) or grand 
unified theories (GUTs) from the heterotic string D □ □ □ □ and superstrings [7j. Partially 
successful results have been obtained. For example, many string theoretical models can explain 



the existence of three chiral families at low energy scales and in principle 

can also provide a natural framework for understanding the fermion masses and mixings ^J. In 
this paper, we construct a new class of three-family models in the heterotic string theory. 

Before presenting our models, it is important to note some caveats common to all known string 
models. They are due to two main difficulties facing the string theoretical model constructions. The 
first difficulty concerns the compactification of the string itself, i.e. the mechanism by which the 
desirable string vacuum is selected. The vacua of string theory compactifications are parameterized 
by many scalar fields with flat potentials. These fields are the modulus fields. They characterize 
the sizes and shapes of the compactified spaces and the strengths of the string interactions; none 
of them can be fixed in perturbation theory [l^. 1 The modulus problem and the related issue of 
supersymmetry breaking will not be dealt with in this article, instead we will simply assume that 
the moduli are fixed by some unknown mechanism at the string scale. 



Stabilizing moduli by fluxes in the context of heterotic string theory has been discussed recently in [T^l . 
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The second difficulty concerns our ignorance of the physics between the electro-weak and uni- 
fication scales. Except for some indication that the SM gauge couplings may unify at about 10 16 
GeV in certain supersymmetric extensions of the SM with minimal matter content jl^], we can 
hardly have any confidence in extrapolating the low energy data by some 14 to 16 orders of mag- 
nitude to the unification scale and infer what the gauge symmetries, matter spectra and physical 
parameters are at that scale. Hence, we are far from having a clear-cut field theoretical model at 
the unification scale to which a string-derived model is supposed to match. Any string construction 
must combine both bottom-up and top-down analyses. 

The new class of string models in this paper are mainly motivated by the recent discussions 
on orbifold GUT models fiol. fsoj] . These GUT models utilize properties of higher-dimensional field 
theories, and have some advantages over conventional 4d GUTs. For example, GUT symmetry 
breaking can be accomplished by an orbifold parity, instead of by a complicated Higgs sector. The 
doublet-triplet splitting problem, which plagues conventional GUTs, can also be solved by assigning 
appropriate orbifold parities to the doublet and triplet Higgs bosons. Note, however, that like all 
field theoretical models in higher dimensions, these GUT models are not renormalizable quantum 
field theories. They can only make sense as low-energy effective theories of some more fundamental 
theory with better ultra-violet (UV) behavior. Our string models provide exactly such kind of UV 
completions, in the sense that they reproduce many interesting features of the orbifold GUTs in 
certain low energy limits. (Connections between orbifold GUTs and an SOio string model have 



already been established in ref. 2l|. The present paper contains more detailed discussions on these 
connections. See also the recent paper [22j.) 

To make the connections between string and field theoretical models more concrete, we consider 
some examples, in particular, the 5d SOio model of ref. and a generalization with bulk gauge 
group Eg. In these models, the extra dimension is taken to be an orbicircle S 1 /^ and the 4d 
effective theory has a Pati-Salam (PS) symmetry, SU4C x SU2L x SU2R The technical app_aratus 
we adopt to build string models is the simplest abelian symmetric orbifold compactification 
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24| of the heterotic string |2(. More specifically, we consider a non-prime-order Zg orbifold 
(or equivalently, Z2 x Z3) model with the orbifold twist vector V6 = g(l,2, —3). To achieve three 
chiral PS families at low energies, we also introduce several (in fact, two) discrete Wilson lines 
[2^1 • 2 It is obvious that the third compactified complex dimension has a Z2 symmetry in the 



Prime-order orbifold models (such as the Z3 orbifold models) with Wilson lines 0, ITol . ITU and non-prime-order 
orbifold models without Wilson lines [2^1 have been extensively studied in the literature. Non-prime-order orbifold 
models with Wilson lines, on the other hand, possess a number of complications, and to our knowledge they 
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model, hence it can consistently be taken to be the root lattice of the SO4 Lie algebra. The string 
models are effectively 5d when the length of one of the SO4 simple roots is large compared to the 
string scale, while all other dimensions are kept comparable to the string scale (i.e. the geometry 
of the compactified space is equivalent to that of the orbifold GUTs, S 1 /Z2). In this limit, the 
heterotic models are similar to the orbifold GUT models in the following respects: 

• The 5d N = 2 supersymmetry 3 is broken to that of N = 1 in 4d by the Z2 orbifold twist and 
the "bulk" gauge group is broken to two different regular subgroups at the two inequivalent 
fixed points by degree-2 non-trivial gauge embedding and Wilson line. The surviving gauge 
group in the 4d effective theory is the intersection of groups at the fixed points. It is the 
PS group in our models. More specifically, we find two types of models. In the first type we 
have an Eg symmetry in the 5d bulk which is broken to SO10 and SUg x SU2 respectively. 
In the second type we have an SO10 x SU2 in the bulk, broken to PS at one of the two fixed 
points. 

• Untwisted-sector and twisted-sector states that are not localized on the Z2 fixed points of the 
SO4 lattice can be identified with the "bulk" states of the orbifold GUT. Interpretation of 
the Kaluza-Klein (KK) towers of the bulk gauge and matter fields agree in the string-based 
and orbifold GUT models. 

• Twisted-sector states that are localized on the Z2 fixed points of the SO4 lattice have no field 
theoretical counterparts, although they can correctly be identified with the "brane" states 
of the orbifold GUT. In the orbifold GUT models, these states are only constrained by the 
requirement of (chiral) anomaly cancellation. 

Of course, string theoretical models are more intricate than the corresponding field theoretic 
orbifold GUT models. They need to satisfy more stringent consistency conditions and thus they 
are physically more constrained. We find it is highly non-trivial (or impossible) to implement all 
the features of the orbifold GUTs. For example, we cannot arbitrarily place the three families of 
quarks and leptons in the bulk or on either brane. Moreover, the very act of obtaining three families, 
along with their respective locations, is fixed by the requirement that the gauge embeddings and 

have not been studied to the same extent. Our work can be regarded as the first serious attempt at constructing 
three-family models from non-prime-order orbifolds. 
3 By N = 2 supersymmetry in 5 or 6d, we mean the minimal number of supersymmetries in these dimensions, (i.e. 
the fermions satisfy the pseudo-reality condition). It reduces to !N = 2 in 4d by dimensional reduction and is 
sometimes called N = 1 supersymmetry in the literature. 
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Wilson lines have to satisfy the modular invariance conditions [J, |26(]. In addition, we cannot 
utilize the orbifold projections to remove all the (6, 1, 1) color-triplet states as in the SOio orbifold 
GUTs j^jj and at the same time obtain three families. We also find many massless states carrying 
unconventional representations under the SM gauge group. These exotic states are commonplace 
in almost all known three-family models. Whether these models can give rise to satisfactory 
phenomenology needs more detailed knowledge of the low-energy effective actions. The present 
status of our analysis is contained in this paper. 

The paper is organized as follows. In sect, [n] we briefly review 5d field theories on the orbicircle 
S 1 /Z2 and present two orbifold GUT models with bulk gauge groups SOio and Eg. The latter 
(model Al 4 ) is a novel 5d model with many nice phenomenological features. Then in sect. Hill we 
discuss the heterotic string construction of model Al. Using this model as a guide we compare 
the heterotic string construction with generic orbifold GUT models by restricting the compactified 
space to a specific type (which is referred to as the orbifold GUT limit). We show the equivalence 
between the matter states (in the untwisted and some twisted sectors) in string-based models and 
the bulk states in orbifold GUTs, as well as their KK excitations. We interpret orbifold parities 
(for the bulk states) in the orbifold GUTs in string theory language, and explain why the gauge 
embeddings and Wilson lines cannot project away all the (6, 1, 1) color-triplet states. These states 
may be needed to break the PS group to that of the SM, as in the field theoretical model of 
sect. Ill Bl In sect. II VI we focus on more of the phenomenological aspects of model Al. In sect. II V Al 
we discuss gauge coupling unification and the determination of the compactification and string 
scales. In sect. UVBl we examine the allowed Yukawa couplings (at both the renormalizable and non- 
renormalizable levels) and their phenomenological consequences, concentrating on the possibility 
of breaking the PS symmetry, mass generation for the color-triplet fields and SM fermions, and 
proton stability. We conclude in sect.0 listing the pros and cons of the present models. Hopefully 
one can learn from the problems to design better models in the future. 

We have made an effort to make the paper more accessible to field theory model builders. Many 
of the details of string constructions are relegated to four appendices. In appendix ^ we review 
the construction of non-prime-order orbifold models with Wilson lines, highlighting its differences 
with the prime-order orbifold construction. In appendix IB II we present three three-family 
models with PS gauge symmetry. The complete matter spectra are listed in appendix IB 21 where 
we also explain the notation for the twisted-sector states. In appendix IC II we review the string 

4 This model is denoted Al since it corresponds to the first of several string models discussed in the paper. 
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selection rules necessary for determining non-trivial Yukawa couplings in a 4d effective theory, and 
in appendix IC 21 list some allowed couplings involving operators of interest in model Al. Finally, in 



appendix 
relations 



I we study gauge coupling unification and derive the Georgi-Quinn- Weinberg (GQW) 
2711 in the orbifold GUT limit. These relations allow us to determine various mass scales 



in our models. 

II. 5D ORBIFOLD GUT MODELS ON S 1 /Z 2 

Let us briefly review the geometric picture of orbifold GUT models compactified on an orbicircle 
S 1 /Z2- The space group of S 1 /^ is composed of two actions, a translation, 7 : x 5 — > x 5 + 2irR, 
and a space reversal, T : x 5 — > — x 5 . There are two (conjugacy) classes of fixed points, x 5 = 2mrR 
and (2n + 1)ttR, where neZ, 

The space group multiplication rules imply 7P7 = IP, so we can replace the translation by 
a composite Z2 action T' = 77 : x 5 — > — x 5 + 2ttR. The orbicircle S 1 /^ is equivalent to an 
K/(Z2 x X' 2 ) orbifold, whose fundamental domain is the interval [0, nR], and the two ends x 5 = 
and x 5 = irR are fixed points of the Z2 and Z 2 actions respectively. 

A generic 5d field <3? has the following transformation properties under the Z2 and Z' 2 orbifoldings 
(the 4d space-time coordinates are suppressed), 

•P : $(x 5 ) -► $(-x 5 ) = P<S>(x 5 ) , y' : <£(x 5 ) ^>(-x 5 + 27ri?) = P'$(x 5 ) , (2.1) 

where P, P' = ± are orbifold parities. In general cases P' 7^ P; this corresponds to the translation T 
being realized non-trivially by a degree-2 Wilson line (i.e., background gauge field). The four com- 
binations of orbifold parities give four types of states, with wavefunctions ^^^(x 5 ) ~ cos(mx 5 /P), 
$ + _(x 5 ) ~ cos[(2m + l)x 5 /2P], $_+(x 5 ) ~ sin[(2m + l)x 5 /2P] and $__(x 5 ) ~ sin[(?n + l)x 5 /P], 
where m G Z. The corresponding KK towers have masses 

m/R for (PP') = (++) , 

(2m + 1)/2P for (PP') = (+-) and (-+) , (2.2) 
(m+l)/P for (PP') = ( — )■ 
Note that only the $++ field possesses a massless zero mode. 

A. An SO10 orbifold GUT 

Consider the 5d orbifold GUT model of ref. j^ . The model has an SO10 symmetry broken to 
the PS gauge group, SU4C x SU2L x SU2R, in 4d, by orbifold parities. The compactification scale 



M c = (7ri?) _1 is assumed to be much less than the cutoff scale. (In string theory the cutoff scale 
is given by the string scale M str i ng .) 

The gauge field is a 5d vector multiplet V = (Am, A, A', a), where Am, o (and their fermionic 
partners A, A') are in the adjoint representation (45) of SOio- This multiplet consists of one 4cl INf = 
1 supersymmetric vector multiplet V = (A^, A) and one 4d chiral multiplet E = ((a + iA^)/\/2, A'). 
We also add a 5d hypermultiplet "K = (<fi, <p c , ip, ip c ) in the 10 representation. It decomposes into 
two 4d chiral multiplets H = (<f), ip) and H° = (4> c ,ip c ) in complex conjugate representations. This 
model has an K = 2 extended supersymmetry. The 5d gravitino ^m = (^mi^m) decomposes 
into two 4d gravitini xjA, rjjz and two dilatini tp^. To be consistent with the 5d supersymmetry 
transformations one can assign positive parities to V>5 — V'i an d negative parities to — V^j 

■ipl + "05 ; this assignment partially breaks N = 2 to N = 1 in 4d. 

The orbifold parities for various states in the vector and hyper multiplets are chosen as follows 
(where we have decomposed all the fields into PS irreducible representations) 



States 


P 


P' 


States 


P 


P' 


y(i5,i,i) 


+ 


+ 


S(15,l,l) 






V(l,3,l) 


+ 


+ 


E(l,3,l) 






V(l,l,3) 


+ 


+ 


S(l,l,3) 






V(6,2,2) 


+ 




E(6,2,2) 




+ 


H(6,l,l) 


+ 




F c (6,l,l) 




+ 


ff(l,2,2) 


+ 


+ 


F c (l,2,2) 







We see the fields supported at the orbifold fixed points x 5 = and ttR have parities P = + and 
P' = + respectively. They form complete representations under the SO10 and PS groups; the 
corresponding fixed points are called SO10 and PS "branes." In a 4d effective theory one would 
integrate out all the massive states, leaving only massless modes of the P = P' = + states. With 
the above choices of orbifold parities, the PS gauge fields and the #(1,2,2) chiral multiplet are 
the only surviving states in 4d. The H(6, 1, 1) and H c (6, 1, 1) color-triplet states are projected 
out, solving the doublet-triplet splitting problem that plagues conventional 4d GUTs. 

B. An E 6 orbifold GUT 

We now consider a novel 5d orbifold GUT with an Eg gauge symmetry. In analogy to the model 
m sect. UlAl we take the 5d gauge field, given by (V, E), in the adjoint representation (78) of Eg. 
In addition to this we add a matter hypermultiplet H(27) + H c (27). 
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We define two orbifold parities 

P = exp(7riQ z /3) x P F , P' = exp[3vri(S - L)/2] x P' F , (2.4) 

which break the Eg via P to SOio and then via P' to PS. Qz is the abelian charge in Eg commuting 
with SOio, normalized such that the 27 decomposes to 16i + 10_2 + l4, and Pp, P' F are appropriate 
discrete flavor charges. (For explicit definition of the parities in the corresponding string model, 
see sect. IIIIBI ) It is easy to obtain the following projections to (++) modes, where the first step 
follows from P alone and the second follows from the subsequent action of P', 

V = 78 -> 45 adjoint of PS , 

X = 78^ 16 + 16^ fl + T, 
27 -> 16 ^ h , 

27^ 10 — /i. (2.5) 

In this equation, we have identified the third family of quarks and leptons as well as the MSSM 
Higgs-doublet pair (h = Hu + Hp where Hjj and Hp are the MSSM Higgs doublets responsible 
for the up- and down-type quark/charged lepton masses), 

/| = (4,1,2), h = (4,2,1), fc= (1,2,2). (2.6) 

As a consequence of the fact that the third family and Higgs doublet come from the bulk gauge 
and 27 hypermultiplets we obtain a gauge- Yukawa unification relation, 

Xt = \ = K = #4d = V&roGUT , (2-7) 

where g id is the 4d gauge coupling constant at the compactification scale. This relation can be 
seen by inspecting the 5d bulk gauge interaction 

Til' 



£ dx 5 ^ 5dJ ff c S^ ^g 4A hflf 3 



where g Ad = g 5d VM^. 

Of course, we then need to spontaneously break PS to the SM via the standard Higgs mechanism. 
This can be accomplished when the "right-handed neutrino" fields in 

X c = (4,l,2), x c = (4,l,2) (2.9) 

obtain non-vanishing vacuum expectation values (vevs) 

(u c ) x c = (V c ) r = M PS . (2.10) 



We already have one such state but we need more (if only for anomaly cancellation) . Consider the 
addition of three more 27 hypermultiplets given by 3 x (27 + 27). Upon applying the orbifold 
parities we find 



3 x (27 + 27) -» 2(16) + 16 + 3(10) -> 2{ X C ) + T + 3(C), 



(2.11) 



where C = (6, 1, 1). We now have a total 2(x c + x c ) fields. Note, with one C, one x c , X c P arr an d 
a superpotential given by 



w = x c x c c + Tx c c, 



(2.12) 



we can give mass to the color triplets and also break PS to the SM along a D- and F-flat direction. 



(The D-flatness condition requires {v c 



(V C ) T and (D 



i/x c 



(D 



then the F-flatness 



condition requires further that one of these vevs, say the second, be zero.) In the end, however, we 
must guarantee that the extra x c , X C i C states obtain mass above the PS breaking scale. 

SO 10 brane SU 6 x SU 2 r brane 



Gauge- V, I (78) 



2 x (16) 



27+ 27 

3 x (27+ 27) 







TlR 



FIG. 1: 5d Eg orbifold GUT model with bulk and brane states. The bulk gauge symmetry is broken to SOio 
on the end of world brane at a; 5 = and to SUg x SU2R at a; 5 = irR. The massless sector of the 4d effective 
theory has a PS gauge symmetry. In addition, the bulk contains four hypermultiplets. and the SO10 brane 
contains two spinor representations, giving rise to the first two matter families. 



But what about the first two families? When constructing an orbifold GUT, one has the 
option of whether to place the first two families in the bulk or on either brane. One of the main 
considerations is to avoid rapid proton decay due to gauge exchange and another is to generate a 
hierarchy of fermion masses. If the compactification scale is much smaller than the GUT scale, say 
M c <C Mqut, then it is not possible to place the first two families on the SO10 brane. It would 
however be fine to place them in the bulk or on the SU6 x SU2R brane, since in the first case the 
families are in irreducible representations with massive KK modes, while in the latter case one 
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family is contained in two irreducible representations (15, 1) + (6, 2), also with massive KK modes. 
In both cases, gauge exchange takes massless quarks and leptons into massive states. Hence there 
is no problem with proton decay. If however M c > Mqtjt then one can place the first two families 
on either brane. Unfortunately, in string theory, we do not get to choose easily where to place the 
families. It is determined by the choice of vacuum. In the heterotic string version of the model 
(model Al in appendix IB")) we find two families sitting on the SOio brane, as in fig. ^ 



III. HETEROTIC STRING CONSTRUCTION OF EFFECTIVE ORBIFOLD GUTS 

In appendix [X] we review the rules for constructing heterotic string models compactified on an 
abelian symmetric orbifold with discrete Wilson lines. Then in appendix [B] we construct three 
three-family orbifold models with two Wilson lines, labelled models Al, A2 and B. We have 
obtained the complete spectra of massless states (plus KK excitations for these models in certain 
limits). As we now show, model Al is the string equivalent to the orbifold GUT in sect. Ill Bl 

The following discussion relies greatly on the notation and discussion in appendices El an d 
FBI Briefly stated, the heterotic string combines a lOd superstring for right movers and a 26d 
bosonic string for left movers. However 16 of the 26 left-moving dimensions are compactified on 
the Eg x Eg root lattice. In order to obtain an effective 4d theory, we compactify six of the remaining 
ten dimensions on a symmetric orbifold defined by a six torus modded by a point group 7h§ with 
the twist vector 

v 6 = ^(1,2,-3), (3.1) 
6 

i.e. the three compactified complex coordinates transform as Z% — > exp(2-7rivg)Zj under the twist. 
The embedding of orbifold twists in the gauge degrees of freedom is realized by gauge twists, V, 
and lattice translations by discrete Wilson lines, W. In abelian orbifolds these vectors simply shift 
the appropriate Es x Es roots. 

To be definite, we choose the six torus as the Lie algebra root lattice G2 © SU3 © SO4, as shown 

in fig. [21 Denoting the basis of the lattice by e = I I © ( 3 I © ( I > wri ° se inner product 

\ e 2 ) \ e 4 ) \e 6 ) 

gives the Cartan matrix of the corresponding Lie algebra, the "Lq discrete symmetry can be realized 
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G 2 SU 3 S0 4 




FIG. 2: Fundamental region of the root lattice G2 © SU3 ® SO4. The filled circles, crosses and squares 
represent fixed points in the T\, T^ t A and T3 twisted sectors. See appendix 151 for further details. 

by the Coxeter element, 5 




under which the basis is transformed to Ce. The Coxeter element has eigenvalues e ±17r / 3 , e ±2l7r / 3 
and e ±17r , thus the three two-dimensional sub-lattices have degree-6, 3, and 2 cyclic symmetries, 
and the corresponding numbers of fixed points are 1, 3 and 4. 

There are three Kahler class moduli (^2,3), whose real parts parameterize the sizes of the three 
tori, and one complex structure modulus (U3), which parameterizes the shape of the third torus. 
Explicitly, ReT3 = 2RR' sin a\ and U3 = § e^, where R, R' are the lengths of the two axes of the 
S04-lattice and cp their relative angle. These moduli are arbitrary parameters. One may make the 
length of one axis (along which one puts the degree-2 Wilson line, W2), say R, large compared 
to the string length scale while keeping all other dimensions small. In this limit (for length scales 
larger than the string scale but smaller than the radius R), the low energy theory is effectively five 
dimensional. 6 The SO4 lattice, on which only the Z2 sub-orbifold twist acts, has four fixed points. 
With only one degree-2 Wilson line, the fixed points split into two inequivalent classes, labelled by 
the winding number ri2 = 0, 1. Thus in our setup the fifth dimension is equivalent to the orbicircle 
S 1 /Z2 where each of the two fixed points has a degree-2 degeneracy. 

Note that we can reinterpret the *Zq models of appendix [B] in terms of the equivalent Z2 x Z3 
orbifold (where the Z2 (Z3) sub-orbifold twist acts on the G2 and SO4 (G2 and SU3) sub-lattices). 
This point of view is more useful for our comparisons with the orbifold GUTs in sect. Ill Bl Labelling 

5 The Coxeter element is an inner automorphism of the lattice, composed of products of Weyl reflections of the 
corresponding root lattice. For example, the Coxeter element of G2 is simply S1S2 where si, S2 are the two 
reflections with respect to planes orthogonal to the two simple roots. A generalized Coxeter element may also 
include outer automorphism of the lattice. 

6 It should be obvious that our construction can be generalized to 6d models, simply by taking both R and R' large 
compared to the string length scale. These models are related to 6d orbifold GUTs compactified on T 2 /Z2. 



11 



a twisted sector in the Zq model by where k = 1, 2, • • • , 5 and in the Z2 x Z3 model by To. a 
where k = 0, 1 and I = 0, 1, 2, then the correspondence between the twisted sectors in the Zq and 
Z2 x Z3 orbifolds is the following: 



Zq orbifold 


T\ 


T 2 


T 3 


T 4 


T 5 


Z 2 x Z 3 orbifold 


T (l,2) 


T (o,i) 


T (l,0) 


2(0,2) 


T (i.i) 



(3.3) 



The T2 t 4 sectors, which will shortly be identified with the bulk states in the language of orbifold 
GUTs, have k = 0, I = 1,2; therefore they are untwisted by the Z2 twist. 



A. Model Al from the Zq orbifold compactification 

We now examine model Al of appendix [Bj Consider first the model with only the Z3 sub- 
orbifolding being imposed (i.e., with twist vector V3 = 2vg, gauge twist V3 = 2Vg and a degree-3 
Wilson line W3, where V6, Vg and W3 are given in eas. l3.HlBll and lB3)l . we find a 6d 3Sf = 2 model 
with observable-sector gauge group Eg (modulo abelian factors). Matter fields of the observable 
sector consist of 6d N = 2 hypermultiplets in the following representations, 

U sectors : 27 + 27, T sectors : 3 x (27 + 27) . (3.4) 

The remaining Z 2 twist acts as a space reversal on the third compactified complex dimension, 
Z3 — > —Z%. The Z3 models have two gravitini with the SOg momentum vectors, r = ^(1, 1, 1, 1) 
and ^(1,-1,-1,1), in the Ramond sector of the right-moving superstring (see appendix lAl for 
notation). Only one of them, r = |(1, 1, 1, 1), satisfies the Z 2 projection, r • v 2 = Z. Hence the 
N = 2 supersymmetry is broken to that of INT = 1 in 4d. 

Gauge symmetry breaking induced by the Z 2 orbifolding is as follows. The twist vector v 2 
is embedded in the gauge degrees of freedom in two different ways, with gauge twists V2 and 
V 2 = V2 + W2 where V2 = 3V@ and W2 is given in eq. IB3I Eq generators in the Cartan-Weyl 
basis are transformed under the Z2 action as £p — > e 27rlP ' V2 £'p and Ep — ► e 2wiP '^Ep, thus the 
linearly-realized gauge groups consist of roots satisfying P • V2 and P • V 2 = Z respectively. The 
pattern of symmetry breaking in the observable sector can be summarized as follows: 

^ SO 10 
E 6 , PS 

SU 6 x SU 2R (3.5) 
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At the final step we have the complete Zg model with two discrete Wilson lines being imposed 
simultaneously; this gives the PS symmetry group in the 4d effective theory. 

In these two inequivalent implementations of the Z 2 twist the non-trivial matter fields of SO 10 
and SU6 x SU2R are: 



Sectors 


so 10 


SU 6 x SU 2R 


Ui 


16 


(15,1) 




10 


(6,2) 


u 3 


16 + T6 


(20,2) 


T (0,l) 


2 x 16+ + 10_ 


2(6,2)+ + (15,l)_ 


T (0,2) 


l6_ + 2 x 10+ 


(6, 2)_ +2(15,1)+ 



(3.6) 



where the subscripts ± represent intrinsic parities, 

V = l4>- (3-7) 

p depends on the twist eigenvalue, 7, and the oscillator phase, tfi; they are defined in appendix lAl 
Note that p = + for gauge and untwisted-sector states, and p = + and — have multiplicities 2 and 
1 respectively for non-oscillator Tm\/T(m) states. 

Massless states in the untwisted and T(o,i)> ^(0,2) twisted sectors of model Al are the intersec- 
tions of those of the SO10 and SUg x SU2R models. This can be seen from the group branching 
rules. For example, the T^i) -sector matter has the following branchings, 

SO10 -> SU 4C x SU 2L x SU 2R 

16+ = (4, 2,1)+ + (4, 1,2)+, 

10_ = (6,1,1)_ + (1,2,2)_, (3.8) 
SU 6 x SU 2R SU 4C x SU 2L x SU 2R 

(6,2)+ = (4, 1,2)+ + (1,2, 2)+, 

(15,1)_ = (4,2,1)_ + (6,1,1)_ + (1,1,1)_. (3.9) 

The states in common, 2(4, 1, 2)+ + (6, 1, 1)_ , agree with that of the T 2 -twisted sector in eq. IB6I 
Massless fields in the other, i.e. TW 2 )(= T%) and Tn \(= T3), twisted sectors are the unions 
of those of the SO10 and SUg x SU 2R models. Therefore there are two sets of states, furnishing 
complete representations of SO10 and SUg x SU 2R respectively. For example, the T\ sector of model 
Al contains (4, 2, 1) + (4, 1, 2) and (4, 1, 1) + (1, 2, 1), they are in the complete representations 
16 of SO10 and (6, 1) of SUg x SU 2R . In the notation of appendix^ these two sets of states have 
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quantum numbers n-2 = and 712 = 1. (These quantum numbers are the winding numbers along 
the direction where the W2 Wilson line is imposed.) The n<i = and n<i = 1 fixed points are thus 
the SO10 and SU6 x SU2R branes in the orbifold GUT language. 



B. Identifying orbifold parities in string theory 



To a certain degree, the above Eg heterotic model gives a string theoretical realization of the 
orbifold GUT in sect. HTBl Better yet, we also achieve an understanding of the orbifold parities in 
terms of string theoretical quantities. Specifically, the analogue of orbifold parities, eq. 12.41 in our 



string models can be defined as follows [21] 

P = p e 27ri(P-V 2 -r-v 2 ) 



P' =;pe 27ri(P-V 2 -r-v 2 ) 



(3.10) 



where V2 and V 2 are the two inequivalent gauge embeddings of the Z2 twist in sect. IIII Al and p 
is the intrinsic parity. 



SUq 




S0 4 
V, £ (78) 



27+27 



3 x (27+27) 

FIG. 3: G 2 ®SU 3 ®S0 4 lattice with Z 3 fixed points. The fields V, E, and 27(e U\) + 27(e U 2 ) arc 
bulk states from the untwisted sectors. On the other hand, 3 x (27 + 27) are "bulk" states located on the 
^(o,i)/^~(o,2) twisted sector (G2, SU3) fixed points. 




kR +2x(1,2)] 



FIG. 4: G2 ® SU3 © SO4 lattice with fixed points. The T^^/T^) twisted sector states sit at these 
fixed points. 



These parities can be deduced from the generalized Gliozzi-Scherk-Olive (GSO) projector 1C, 
28(, as in the paragraphs after eq. IA141 Since the terms in the exponents, P • V2 — r • V2 and 
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G 2 SU 3 S0 4 

SO 10 SU 6 xSU 2R 




71R 



FIG. 5: G2 © SU3 © SO4 lattice with Z 2 fixed points. The T(i : o) twisted sector states sit at these fixed 
points. 

P • V 2 — r • V2, take integral or half-integral values, P and P' are either + or — . The orbifold 
translation corresponds to the difference in P and P', i.e. T = e 27rlP W 2 The p ; p 1 an j j 1 [ n string 
models have exactly the same properties as that of the orbifold GUTs. 

Evidently, in the E 6 orbifold GUT model of sect. Ill Bl states supported at the SO10 and SUg x 
SU2R branes are those with parities P = + and P' = +, and states in the 4d effective theory 
are those with parities P = P' = +; this agrees with the string theoretical interpretation, since 
the parities in eq. 13.101 are nothing but the required GSO projections for the gauge, untwisted 
and T(oi)/T(o2) sector states (i.e. the bulk states) in string models. (The massless states, i.e. 
P = P' = + modes from bulk and JV^/Im^) twisted sectors are shown in fig. El) From 
information gathered in sect. MI Al and appendix El we can also deduce the P and P' parities for 
the various bulk matter states. They are listed in table [IJ 

KK masses for these bulk states can also be derived in string models. The mode expansions of 
the coordinates corresponding to the SO4 lattice are X l L R = x l L R + p l L r {t ± a) + oscillator terms, 
with p l L , p R given by eq. IA3I The Z2 action maps m to — m, n to — n and W2 to — W2, so physical 
states must contain linear combinations, \m, n) db | — m, —n); the eigenvalues ±1 correspond to the 
first Z2 parity of the orbifold GUT models. The second embedding corresponds to a non-trivial 
Wilson line; it shifts the KK level by m — > m + P • W2. Since 2W2 is a vector of the integral 
Eg x Eg lattice, the shift P • W2 must be an integer or half-integer. In the orbifold GUT limit when 
the winding modes and the KK modes in the short direction of SO4 decouple, eq. IA.3I reproduces 
the field theoretical mass formula in eq. 12.21 

As seen in sect. IIII Al matter states in the Tn^/Tnp) an d ^(1,0) twisted sectors, which may 
be identified with the first two families, are localized on the two inequivalent fixed points in the 
SO4 lattice. They are the SO10 and S\Jq x SU2R brane states (See figs. UlandEl). These twisted- 
sector states are more tightly constrained than their orbifold GUT counterparts. In orbifold GUT 
models the only consistency requirement is the chiral anomaly cancellation, thus one can add 
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TABLE I: Parities for the bulk states in model Al, computed from eg. 13. 101 The states have been decomposed 
to the PS irreducible representations. 



Multiplicities 


States 


P 


P 1 


States 


P 


P' 


1 


7(15,1,1) 


+ 


+ 


£(15,1,1) 




~ 


1 


7(1,3,1) 


+ 


+ 


£(1,3,1) 




— 


1 


7(1,1,3) 


+ 


+ 


£(1, 1,3) 






1 


7(6,2,2) 


+ 




£(6,2,2) 




+ 


1 


7(4,2,1) 




+ 


£(4,2,1) 


+ 




1 


7(4,1,2) 






£(4,1,2) 


+ 


+ 


1 


7(4,2,1) 




+ 


£(4,2,1) 


+ 




1 


7(4, 1,2) 






£(4,1,2) 


+ 


+ 


1 


#(4,2, 1) 


+ 


+ 


# c (4,2,l) 






1 


#(4,1,2) 


+ 




# c (4,l,2) 




+ 


1 


#(6,1,1) 




+ 


# c (6,l,l) 


+ 




1 


#(1,2,2) 






# c (l,2,2) 


+ 


+ 


2 


#(4,2,1)+ 


+ 


— 


# c (4,2,l) + 


— 


+ 


2 


H(4, 1,2)+ 


+ 


+ 


# c (4,l,2)+ 






2 


#(6,1,1)+ 






# c (6,l,l)+ 


+ 


+ 


2 


#(1,2,2)+ 




+ 


# c (l,2,2)+ 


+ 




1 


#(4,2, 1)_ 




+ 


H c (4,2,l)_ 


+ 




1 


#(4, 1,2)_ 






# c (4,l,2)_ 


+ 


+ 


1 


#(6,1,1)- 


+ 


+ 


# c (6,l,l)_ 






1 


#(1,2,2)_ 


+ 




ff c (l,2,2)_ 




+ 



arbitrary numbers of vector- like representations to the branes. String models have to satisfy more 



stringent modular invariance conditions |J, |2(| (of course, one-loop modular invariance guarantees 
the model is anomaly free, up to a possible anomalous abelian factor 2[|), which also constrains 
any additional matter in vector-like representations. 



Other models 



In this subsection, we discuss the other two models in appendix[B]in the orbifold GUT language, 
for completeness. These models do not have matter-Higgs couplings at the renormalizable level, 
and may have limited phenomenological interest. 

Model A2 has already been analyzed in ref. |2j|. In the 5d bulk, it has an SOio x SU2 gauge 
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symmetry and the following set of matter states, 

[/sectors : (16, 1) + (1,2), T sectors : 3(16, 1) + 6(10, 1) + 15(1, 2) . (3.11) 

The bulk gauge group is unbroken at the fixed point rr 5 = and broken to the PS group at 
x 5 = irR respectively, the states supported at these two points are 



Sectors 


SOio x SU 2 


PS 


Ui 


(16,1) 


(4,2,1) 


U 2 


(1,2) 


(4,1,2) 


lh 




(6,2,2) 


T (o,i) 


2(16,1)+ + 2(10, 1)_ 


(6,1, 1)_ +2(1,2,2)+ 




+2(1,2)+ + 4(1,2)_ 


+2(6, 1,1)+ + (1,2, 2)- 






+(4,2,1)_ + 2(4,1,2)+ 


T (0,2) 


(16, 1)_+ 4(10,1)+ 


2(6, 1,1)+ + (1,2, 2)_ 




+8(1,2)+ + (1,2)_ 


+(6, 1,1)_ +2(1,2,2)+ 






+(4,1,2)_+2(4,2,1)+ 



(3.12) 



Model B is similar to model Al, with an Eg bulk group and the same set of bulk states as in 
eq. 13.41 The Eg group is broken to SOio and SUg x SU2L respectively at the two fixed points, and 
the matter states are 



Sectors 


SOio 


SUg x SU 2L 


Ui 


16 


(6,2) 


lh 


10 


(15,1) 


u 3 


16 + 16 


(20,2) 


T (o,i) 


2 x 16+ + 10_ 


(6, 2)_ + 2(15,1)+ 


T (0,2) 


16„ + 2 x 10+ 


2(6, 2)+ + (15,1)- 



(3.13) 



The 4d effective theories of models A2 and B have a PS symmetry, and the complete matter 
content is listed in appendix El Similar to model Al, matter fields in the untwisted and T2/T4 
twisted sectors can be traced back to the states in the above two tables, by using appropriate 
group branching rules. We also note that both models contain two families of chiral matter from 
the fixed point at x 5 = and one family from the bulk, a common feature to all our models. This 
feature predicts a non-abelian dihedral D4 family symmetry - a novelty in string model building - 
as we will see in the next section. 
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D. The color-triplet problem 



A major motivation for constructing orbifold GUT models is the well-known doublet-triplet 
splitting problem in conventional 4d GUTs. Orbifold GUTs solve this problem by assigning appro- 
priate orbifold parities to the Higgs doublets and triplets, such that the triplets are automatically 
projected out of the effective theory Q|. 7 We have already seen this in the SOio model in sect. Ill Al 
This mechanism, however, usually cannot be trivially implemented in heterotic models. 

The difficulty is largely due to the intricate nature of string models. These models need to satisfy 
delicate modular invariance consistency conditions P, [3| an d are physically more constrained than 
the orbifold GUTs. Before imposing any Wilson line, the Zg models of eqs. IB II and IB2I always 
contain the 10 representation of SOio simultaneously in several sectors. We find it is impossible 
to design modular- invariant Wilson lines to fulfill the following requirements: (a) break the gauge 
group to PS in 4d, (b) give rise to three chiral families, and (c) eliminate the color triplets altogether. 
Furthermore, the 10 representation of the T3 sector in model Al does not suffer from additional 
projections even when the W2 Wilson line is turned on. Indeed, it simply decomposes to the 
(6, 1, 1) + (1, 2, 2) representations under the PS group. 

Although the presence of many color triplets is a nuisance, one 3 + 3 pair may be necessary to 
facilitate the breaking of PS to the SM gauge group, as illustrated in the Eg model in sect. Ill Bl 
Moreover, it is not entirely clear that the color triplets in our models pose the same problems 
as in conventional GUTs. Indeed, although there are color triplets (those of the T1/T3 twisted 
sectors) with doublet companions having exactly the same quantum numbers, in general we also 
have (6, 1, 1) and (1, 2, 2) states with different quantum numbers, in all three models. The usual 
doublet-triplet problem does not necessarily apply for the second situation. We need to check 
whether it is possible to make all color triplets sufficiently heavy while (hopefully) keeping one 
MSSM Higgs-doublet pair light; this requires a better understanding of the effective actions of our 
models and will be examined for model Al in the next section. 

IV. PHENOMENOLOGY OF MODEL Al 

We have seen in the previous sections that the Eg orbifold GUT and heterotic Al model match 
nicely in the low energy regime. In this section, we study some phenomenological issues for the Al 
model. We first study gauge coupling unification, relying on a simplification due to the correspon- 

7 Of course, more conventional field theoretical mechanisms [3pl | have been widely studied in the literature. 
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dence between the field and string theoretical models. We then study Yukawa couplings (including 
both renormalizable and non-renormalizable couplings), concentrating on several immediate phe- 
nomenological questions: (i) breaking of the PS symmetry to that of the SM, (ii) mass generation 
for the color triplets, (iii) proton stability, and (iv) matter-Higgs Yukawa couplings. These cou- 
plings are introduced by hand in the orbifold GUTs. In string models we no longer enjoy the 
same kind of freedom. In fact, these couplings are determined by string selection rules (reviewed 
in appendix IC The low energy phenomenology also depends crucially on the flat directions of 
the effective N = 1 model. However, we do not attempt to solve this complicated problem here. 



A. Gauge coupling unification and proton decay 

As discussed earlier, since the first two families are located on the SOio brane, proton decay 
constraints require that the 5d compactification scale M c be greater than O(10 16 ) GeV. However 
all these GUT scale thresholds must be consistent with low energy gauge coupling unification. 
Consider the solution to the 5d renormalization group (RG) equations, i.e., the GQW equations, 
in the orbifold GUT limit, 8 given by 

2?r - ■ 2?r ■ OlSSM i „ M PS . /, PS , j, n i „ Mitring 



+ bg ^ + (b£ + 6 brane ), log ■ 



Oii(n) ttstring A* Mp S 

where Mps is the PS breaking scale and a s tring is the gauge coupling at the string scale. In addition, 
in the weakly coupled heterotic string we have the boundary condition 

_ 7T / xvLpi ) 2 

^string 

where the first term is the tree level result and the second is a universal one loop st ring y correction. 



EfJ^y + iA"* (4.2) 
4 \M strin J 2 y ' 



361 ] we see that 



The latter correction depends on the value of the T3, U3 moduli. Following ref. 
A umv is a finite function of its argument (with a mild singularity when T3 = U3, modulo PSL(2, Z) 
transformations). Since the universal correction is not significant, we use the tree level formula in 
the following. 



In principle these equations can be derived from a string theory calculation, following refs. 

EJ El IH Hi 0. 

However, it is difficult to obtain the GQW equations in analytic form in string models with discrete Wilson lines 
[35L which makes the calculation less practical for our purposes. Instead, in deriving eg. 14.11 we have worked in 
the orbifold GUT limit, and assumed the most important contributions to the gauge threshold corrections come 
from the KK tower of the large dimension of the SO4 lattice, with a physical cutoff at the string scale, M s trin g . 
See appendix |H| for more details. 
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Eq. 14.11 can be compared mathematically to the 4d equations given by 

271 - ^ +&r SM log^+6fe, (4.3) 



Qti(ll) ^GUT ^ 
where Mqut — 3 x 10 16 GeV, a GUT — ^4 an< ^ we have included a threshold correction at MquTj 
required in order to fit the low energy data. 

With the bulk field parities given in table HJ we find the beta function coefficients, ^(fr++ + 
b^_) = (1, 1, 1), = (|, —3, —3). The brane contributions include that of the two PS families, 

frbrane = (4,4,4), and those from extra matter fields, n' Q 1, l) + mo(l, 1, 1) + n 2 fl (ljj> 0> 0) (with 
n' 6+s < 4, nio < 6 and n' 2R < 12). Equating the difference 2-k/olz(h) — lixja^ip)^ eqs. 14. II and 14.31 
gives 

M PS ~ e- 3/2 M GU T ^ 7 x 10 15 GeV. (4.4) 

Then equating the difference 27r/a2(/u) — 27r/ai(/i), we find 

/Vf . . 8 — 3n' 9 „ + 6n' . 

log ^IHil ~ ^ «±e_ (4.5) 

S M GUT 32 + 2n^-4n' 6+5 ' 1 7 

which results in a maximum value for M s t r ing for n' 2R = 0, ?ig + g = 4, given by 

M2^ g oi e 2 M GUT ^ 2 x 10 17 GeV. (4.6) 

Finally we have 



"gut _ i ^ a GVT 

^string 27T 



M GUT , / M string . 3 

lo g — T? ( n io + n„ . «) log — & + - 



(4.7) 



M c y 6 + 6; ^ * Mqut 2 

Using M str i n g = M s ™^ g and eq. 14.21 for the tree level heterotic string boundary condition, we find 
there is no solution, consistent with M s t r i n g > M c ~ M G ut > Mps- The problem is that the 
value of a s tring gi yen by ea. 14.21 is much too small (a s tring "C ocut) and it cannot be obtained by 
logarithmic running above the compactification scale (note, 6 Ee = and thus there is no power- 
law running). This problem suggests that non-trivial (perhaps non-perturbative) string boundary 
conditions are required for consistency. 

We have considered the lid Hofava-Witten extension 
boundary condition given by j^] 



371 ] of the perturbative heterotic string 



2vr 1 /M P1 X 2 



(4.8) 



^string 2(4vr) 5 /3M, / 9 V M * 
where M* is given in terms of the lid Newton's constant by k 2//3 = M~ 3 and p is the size of the 
eleventh dimension. Now using ea. 14.81 we find solutions for M s t r ing — M* = 2M G ut> M c ~ Mps ~ 
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e _3//2 MGUT with n' 2R = n ' 6+ Q = 4 and M*p ~ 2. Of course, this solution provides an enhanced 
proton decay rate due to dimension-6 operators with the dominant decay mode p — > e + 7r°. The 
decay rate for dimension-6 operators is given by j^] 



Afv W0.015GeV 3 \ 2 



r( p^V),i a >ir| „ GeV l^w^l 



73- 2 



3xl0 33 , 0.015 GeV ^ 



/^lattice 

where /lattice i s an input from lattice calculations of the three quark matrix element. 9 Recent 
results give a range of central values /lattice = 0.007 — 0.015 |40J. Note, the present experimental 
bound for this decay mode from Super-Kamiokande is 5.7 x 10 33 years at 90% confidence levels 
4l| . Thus this prediction is not yet excluded by the data, but it should be observed soon. 



B. Yukawa couplings 

1. PS symmetry breaking, mass generation for color-triplets and proton stability 

To successfully break the PS symmetry to that of the SM and generate mass for unwanted 
color triplet states, the 5d Eg heterotic model should contain non-trivial couplings of the form in 
eq. 12.121 The model, however, contains additional color triplets. They could, in principle, develop 
mass through non-trivial Yukawa couplings to, say, singlet fields. In order to verify if this is a 
possibility, we need to know whether the required couplings exist in the 4d effective theory of the 
string model. For this purpose we are particularly interested in non-trivial couplings containing PS 
invariant operators, (6, 1, 1)(6, 1, 1), (6, 1, 1)(4, 1, 2)(4, 1, 2) and (6, 1, 1)(4, 1, 2)(4, 1, 2), that are 
allowed by string selection rules. (In field theory, all gauge invariant operators would be allowed.) 
These rules are reviewed in appendix IC H and the relevant operators are given in eqs. IC13I - IC16I 

Cubic, renormalizable, couplings in model Al are determined in eq. IC111 they contain the 
following operators of interest (we label the fields according to table Hi]) . 

W< 3 > d S 1 (C 3 ) Aa (C 3 ) Bp + + (C 3 ) Aa f c BX p + maX c pfi + (C,) a f c A f B , (4.10) 

where a,f3 = 1,2 labels the two 7 = 1 eigenstates in the T2 ; 4 sectors, A,B = 1,2 indicate degen- 
eracies associated with the n' 2 winding number (which corresponds to a hidden S2 permutation 
symmetry). The string selection rules require a + (3 = mod 2, A + B = mod 2, and for the last 

9 To obtain this result we have taken the decay rate for SU5 pffll ) and multiplied the amplitude by an additional 
factor of two to account for the extra gauge exchange present in SOio- 
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term a is unrestricted. Apparently these couplings are not sufficient to break the PS symmetry 
and give mass to all the color triplets contained in the (6, 1, 1) and (4, 1, 2) + (4, 1, 2) states. Thus, 
to achieve our goal, we must also take into consideration higher-dimensional operators. Some of 
them are listed in appendix IC 21 (Of course, there are many more operators with even higher 
dimensions. It is not obvious that "stringy zeroes" exist. We omitted these operators. They may 
or may not disrupt the following discussion.) 

A few observations are appropriate. The fields Cj, i = 1, • • • ,4 transform as a 3_ 2 / 3 + 3 2 /3 
under SU3C X Uiy or as D° + D where D has the quantum number of an anti-down quark. These 
color triplets have multiplicities [in brackets] , Ci [1] , C^p], Ca[4], C?4[2]. In addition C 2 , C3 appear 
in complete SO10 10-plets. We also have states contained in Xai X12 with quantum numbers of 
D, D c , U, JJ C , E, E°. Finally we have the color triplet states in q\, q 2 , ~q~i, q 2 with multiplicity 
2 each. These latter are exotic states with fractional charge ±1/2 for the color singlets and ±1/6 
for the color triplets. 

For the exotic states we find operators of the form q 1 (q\ + q 2 ) multiplied by products of singlets 
S n (up to sixth order), but no operators of the form q 2 {q\ + Q.2)S n to order n = 9. Hence q 2 and 
one linear combination of q\, q 2 remain massless. This is a serious problem for the model, since 
these states are absolutely stable and should have been observed. It remains to be seen whether 
the operator q 2 (qi ± Q2)S n is generated at order n > 10 or if it is forbidden by the string selection 
rules to all orders. 

Now consider the fields Ci, i = 1, • • • , 4 and x%i Xi %■ ^ n this sector we need to both find a way 
of spontaneously breaking PS to the SM, as well as giving all color triplets mass. A related issue is 
the potential problem of rapid proton decay mediated by these color triplets. In particular we must 
eliminate or greatly suppress the following baryon/lepton- number violating effective operators 

fff c (x c S n ) =► QLD + LLE, f c f c f c ( X c S n ) =► UDD. (4.11) 

Note we have checked that these operators are not generated prior to integrating out the color 
triplets, for n < 3. 10 Nevertheless there is a danger that they will be generated in the effective 
theory below the color triplet mass. In fact, consider the renormalizable couplings in eq. 14.101 It 
is evident that an effective mass term of the form (5 , )(C , 4) 2 combined with the coupling C^Xafs + 

10 Of course, it would be better check to any order in n, or better yet, find a symmetry which forbids them to all 
orders. 
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f%f%) leads to the effective operator of the form 

ToT XafiflfB = 





(4.12) 



Similarly, additional baryon/lepton- number violating operators are obtained from an effective C3C4 
mass term. These baryon-number violating operator may be phenomenologically acceptable if the 
coefficient is sufficiently small. However it seems prudent to eliminate the offending mass terms 
by choosing a vacuum configuration where the appropriate scalar vevs vanish. For example, given 
the superpotential terms in appendix C, eq. IC13| we demand that the following vevs (i.e. the 



coefficients of (C4) 2 and C 3 C 4 ) vanish S 2 S 24 + §f ^ = + S2S9S24 + S i § ? ] + SigS 25 §f ] = °- 



In the following we consider the possibility of obtaining a baryon/lepton- number conserving low 
energy effective theory. 11 As a possible proof of existence, we suggest an ansatz where the only 
singlets with non-vanishing vevs are 



where the curly braces represent classes of singlets with the same transformation properties under 
all string symmetries and we use square brackets when we explicitly present a finite set of fields in 
the same class. The corresponding superpotential from appendix IC 21 is 



W D [{S 1 }(C 3 ) 2 + {S 2 S 22 S { v ) }C 2 C 3 + {S l0 sf ) }C l C A j 
+ (ctfffi + c 3 r AX c a + C 4 ( X C J3 + flffi) 

+ (^{S 2 S 24 S 2 Q}C 4 (xi) 2 + [SioS^ + SqSjS^Sis + SiSioSj^ + ^si^JC^xi) 2 

+ [S 2 S 3 S 12 § { 3 2) + S 7 S 14 S 16 (S 6 S 2e + s 2 s 27 )}Ctxlxi) + {S 2& }xVxl- (4.14) 



Note the coefficient of the C 2 C 3 term and the first term linear in C 4 vanishes due to the vevs we 
have chosen, but there may be other higher-dimensional terms which replace them. For example, 
the first element in the second term linear in C4 also vanishes, but the other terms may be non-zero. 
The following vevs are assumed to vanish. 



The conventional wisdom of field theory is to use an R-parity (or family reflection symmetry) to eliminate these 
baryon/lepton- number violating operators. (The R-parity has a bonus of predicting a generic stable neutral 
fermionic superpartner, which makes it even more appealing phenomenologically.) Although from the start our 
string models contain several discrete R symmetries at the level of 4d effective action, it is not clear a priori whether 
any of them can survive symmetry breaking. Note that in our models an unbroken R parity (which is capable of 
distinguishing the Higgs fields, \ c \ from the matter, f c ) does not exist because both 6*4X0/3 and CifXf% couplings 
are allowed by string selection rules at the renormalizable level. 



{S\}, {S10S4 }, [5 , 6-S , 7S'i45i8], [SrSi4Sie(SeS 2 6-\- S 2 S 2 7)], {S 2 q} 



(4.13) 



[S 2 S M + §Y^2] 



[§f ] + S 2 S 9 S 24 + SiS { 2 2) + S i9 S 25 § {2) ] 
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{S2S9S22}, [Sl0&3 + S^S^ + 5'l'S'lo5 , 12'S'32 + SgSxoS^ ], 

{Sg} = Sg + S1S10S21S22 + -Sfl + - - - j [Sg^ + 'SitKSklScj 2 '' + Si3§^)], 

[S 2 2 ^ + 526§ 2 3 ^]' [^S^ + 5 , i3(S , io§3 2 ' ) + Si3§ 2 ^)], 

{5'iO'S'i3'S'2i'S'22}, {5 , gS , io5'i25 , 32}, {S'loSi' 3 '*} (4-15) 

With this choice we guarantee, at least to the order we have checked, that we do not generate 
baryon/lepton- number violating operators in the low energy theory obtained by integrating out 
the color triplets. A self-consistent solution to the necessary set of vevs is given by 

S4 = Sg = Sn = S13 = S17 = S21 = 5*24 = 5*25 = S29 = = S32 = 0, (4-16) 

and all other vevs non-zero. In addition we require [S5/S33 + SicS^] = 0, which may or may not 
require fine-tuning. Unfortunately we are not able to identify a symmetry which would extend this 
result to all orders in string perturbation. This is a serious problem for the model. 

The first term in parentheses of eq. 14.141 gives mass to the color triplets C2, C3, C\ and one 
triplet in C4. Since the doublets /i 2 and /13 have exactly the same quantum numbers as that of 
C2 and C3, they acquire the same mass and also decouple from the low energy effective theory. In 
this way, we obtain just one doublet field h\ (from the C/ 2 sector). It is a good candidate for the 
MSSM Higgs-doublet pair. 

The last term of eq. 14.141 is in the form of eq. 12.121 with two pairs of x° + X c fields. With 
(526) / an F-flat solution is found with (xl) = (xf) = 0, then this part of the superpotential 
reduces to ea. 12.121 Employing F and D flatness conditions, the potential has a flat direction along 
the right-handed-neutrino direction of \\ + X2 ' as i n e( l- I2.1UI The non- vanishing vevs break the 
PS to the SM gauge group and subsequently gives mass to the remaining massless color triplets in 
C4 and the D c + D° components in \\ + X2 ■ The remaining charged states in x\ + X2 obtain mass 
via the super-Higgs mechanism. Lastly, the x\ + Xi a l so acquire mass via the (S26) vev 111 e q- 14.141 

The remaining central problem is whether the required singlets, such as those in eq. I4.13( could 
develop the appropriate vevs along F- and D-flat directions. There are five abelian factors in 
the Al model, one of them is anomalous (which is cancelled by the generalized Green-Schwarz 



mechanism 
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, 4j|, as usual). In general, the anomalous Ui factor destabilizes the original vacua 
and contributes vevs to some of the singlet fields. It requires further investigation to determine 
whether our assumptions in the previous paragraph are substantiated, and moreover whether all 

the abelian symmetries may be broken with the vacuum solutiom (The K = 1 supersymmetry, 

however, is generically preserved in the effective theory.) Refs. have already obtained some 
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necessary conditions for analyzing non-trivial singlet vevs. We leave a detailed investigation for 
the future. 



2. D4 family symmetry 

Before discussing the Yukawa matrices for quarks and leptons, we consider the family symmetry 
of model Al. The third family is a bulk field, while the first two families are located on the two 
Z2 fixed points in the SO4 torus with an SO10 gauge symmetry. One family sits at each fixed 
point (see fig. H}. Since the Wilson line in the SO4 torus lies in the orthogonal direction to these 
two fixed points, the theory is invariant under the permutation of the first two families, labelled 
by an index n' 2 = 0, 1 (or A = 1,2). In addition, the string selection rule, eq. IC81 requires that 
every effective fermion mass operator include an even number of fields with n 2 = 1. Hence these 
effective operators are invariant under a Z2 parity n' 2 — ► —n' 2 . The two operations are generated 

1 

by the two Pauli matrices o\ = I I and 03 = I I acting on a real two dimensional 




-1 

vector. The complete set of operations closes on the discrete non-abelian family symmetry group 
D4 = {±1, ±<7i, ±03, Tio^}- Note that the eight-element finite (dihedral) group D4 is the 
symmetry group of a square. It has five conjugacy classes and five faithful representations. The 
character table is 



Classes 


I 


—I 


±<T1 


±£7 3 


Ti<72 


Doublet - D 


2 


-2 











Singlet - Ai 


1 


1 


1 


1 


1 


Singlet - B 1 


1 


1 


1 


-1 


-1 


Singlet - B 2 


1 


1 


-1 


1 


-1 


Singlet - A 2 


1 


1 


-1 


-1 


1 



(4.17) 



In our models, the first two families transform as the doublet, while the third family transforms as 
the trivial singlet. 

We have many SO10 singlets in our models, transforming as doublets under D4. They appear 
in effective higher dimension fermion mass operators. Consider, for example, two doublets under 
D4 given by {Sa,Sa}- Then in terms of these two doublets we can define bilinear combinations 
transforming as {Ai, A 2 , B±, B 2 }. We have 

iSi<Si + S2S2 ~ A\ 
S1S2 — S2S1 ~ A2 
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S±S 2 + S 2 S± ~ B\ 

SxSx - S 2 S 2 ~ 5 2 (4.18) 

The effective Yukawa couplings are then constructed in terms of D4 invariants. Define the D4 
doublet left-handed quarks and leptons (4, 2, 1) [= f A ] and left-handed anti-quarks and anti- 
leptons (4, 1,2) [= f A \ for the first two families and the Higgs multiplet (1,2,2) [= h]. We then 
have the PS and D4 invariants: 

hA 1 (f 1 f^ + f 2 fi) = hA 1 (f A f%) 

hA 2 (fifZ-f 2 ft) 

^i(A/ 2 c + / 2 /i c ) 

hBzihfl- f 2 f 2 ) (4.19) 

We can also have operators of the form 

h(f A S A )(f B S B ) = M/i/i^i + hflS\ + (/x/ 2 c + hfl)SiS 2 \ (4.20) 

Unfortunately there are, in principle, several possible ways of constructing D4 invariants. We are 
not able to determine, without further string calculations, how to contract the D4 indices. In the 
following we assume, for illustrative purposes, that only the simplest invariants, A\ and B±, appear 
in the effective Yukawa couplings. 

3. Fermion masses 

The only Yukawa coupling in model AI present at leading order is for the third family, given 
by the first term in eq. IC11I From discussions in sect. Ill 131 we conclude that this coupling unifies 
with the GUT gauge coupling at the 5d compactification scale, as in eq. 12.71 Yukawa couplings 
for the first two families come from higher-dimensional non-renormalizable operators. In addition 
they are constrained by the D4 family symmetry. In principle, there are at least two possible types 
of operators. The first type involves operators of the form hi(f 3 f A , f A f§, f A f B ): multiplied by 
suitable singlets, and the second type also involves composite singlets x c X°- The second type of 
operators is particularly important since it has the potential to discriminate up-type quarks and 
charged leptons from the down-type; this is necessary to obtain a realistic CKM matrix and also 
resolve the "bad" GUT relation m s /nid = m^/m e . 

We define the following two composite operators 

Oi=riXa, ®2 = %Xl, (4.21) 
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where the group indices are arranged in all possible ways. From the string selection rules, it is 
straightforward to show that the Yukawa matrix is (we only keep representative terms, see eq. IC22I 
for more complete expressions), 



/ in c(3,9,12) . c (10,22,22,23) ,,, £,(3,9,12) . c (10,22,22,23) ,„ ,„ c (3,12) c (9,10,22,22,23) 



(fif2h)hi 



where 



O2SH 



(3,9,12) 



-,(10,26) 



(10,22,22,23) 



o 2 s [ c 



(3,9,12) Q (10,22,22,23) 



5^ 



+ S, 

(10,26) 



Oi0 2 ^ 3,12> + s ( 



(9,10,22,22,23) 



/I 



(4.22) 



5 W>,-) = J- S a A S b B 

Y^, A=even 



g(fifi,..) = J- S a A S b B ---, (4.23) 
£ A=odd 

with are the family indices of the corresponding singlets (of the sectors). Several comments 
on eq. 14.221 are now in order. 



The structure of the Yukawa matrix is determined by a D4 family symmetry. (See the caveat 
at the end of Section HVEH ) 



• Given the superpotential for color triplets, eq. 14.141 an F-flat direction requires (xl) = 
which gives Oi = 0. If however there is a higher-dimensional operator of the form 5 , 26 , S'x»X2' 
then the combined terms S^XaXi + SxaX.2) has an F-flat solution with (xi), 7^ ® an d 
thus Oi, O2 7^ 0. We will analyze the more general case. 

• Given the superpotential for color triplet masses, our previous solution eq. 14.161 requires 
(Sg.) = 0. Hence the composite operators S^ 3 ' 9 ' 12 \ 5(940,22,22,23) yajjjgh However, it is again 
possible that these terms may still be present when higher order products of operators are 
considered. 

• It is crucial to understand how the PS group indices are contracted and what the cor- 
responding Clebsch-Gordon (CG) coefficients are. In orbifold models, the massless matter 
fields correspond to the (integral) highest weight representations of the level-one Kac-Moody 
algebra. In principle one may extract the desired information from the conformal blocks. 
We have not attempted to perform such a string theoretical analysis. Instead we shall adopt 
a simpler field theoretical approach, following ref. j^ . 

Our aim is to examine phenomenological implications of eq. 14.221 in a simple setting. We shall 
consider two simple cases in the following. 

Case A - 
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First neglect the (13) and (23) entries of eq. 14.221 (it may be reasonable to do so, because they 
are higher order terms), and consider the 2x2 sub-matrix corresponding to the second and third 
families, i.e. 



£(3,9,12) + ^(10,22,22,23) Q 
£(10,26) 1 



(4.24) 



We may take O2 to be in the form of the W operator of ref. [45|. This operator has a vanishing 
(non-vanishing) CG coefficient for the up (down) type fields. One may require 

<0^ 3 > 9 ' 12 )) „ ^± „ A 2 , (5(10,22,22,23)) ^ ^ A 3 ; (425) 

m b m t 

where A ~ 0.22 is the Cabibbo angle. The mixing angle V c b is then approximately \ 2 (So' 26 ^}, 
implying 5^ 0,26) ~ 0(1). 

Next consider the (2 x 2) sub-matrix corresponding to the first and second families. Note that 
this part always has the following form, 



(4.26) 



The democratic form may lead to realistic values for quark masses and mixings. Taking a U) d = 

(3 9 12) 

&u,d(l + £u,d) with Ed « A (which implies an approximate Z2 symmetry between vevs, 02<S'e ' ' = 

m 0(3,9,12)/-, . \ 0(10,22,22,23) r,(10,22,22,23)/-, . \ , c (10,26) c (10,26)x ■ 

02^o (l + £d))"e — 60 (1 + £u) and be —bo ), we obtain the mass 

ratio md/m s and CKM angle V us at correct orders. More suppressed value for 10" d , 
is required for m u /m c . Finally, it is also possible to obtain correct mass relations for the charged 
leptons, rrifi/mr ~ m s /m& and m e /m M ~ md/m s . 

Case B - 

Consider now the case that the (13) and (23) entries are not negligible. Let us parameterize 
the (23) entry in the down sector by and assume the corresponding entry in the up sector is 
smaller (or comparable). For simplicity, we also assume 

5(10,26) = Q The ^ x 2 j sub _ matrix of 

the second and third families is 



2 sf 9 ' 12) + si im23) & 



(4.27) 



We can obtain correct mass ratios m c /mt and m s /mb and mixing angle V c b if J/23 ~ 02<Sc 3,9 ' 12 ^ ~ A 2 

(10 22 22 23) o 

and be ' ' ' ~ A . The discussion on the first and second families follows essentially in the 
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same way as in case A. However, we now need to tune C^S^ 3 ' 9 ' 12 ^ + S^ 10 ' 22 ' 22 ' 23 ^ appropriately. 

Finally consider neutrino masses. The Dirac neutrino mass matrix has the same form as that 
of eq. IC22I Effective Majorana neutrino masses are obtained in eq. IC241 where the non-trivial 
effective operators have the form fafbX^Xj ( a > b = 1, 2, 3, i,j = 1, 2) with suitable powers of singlets. 
The non-vanishing vevs of xl project out the right-handed neutrinos in /|, One then obtains 
a Majorana mass of order < M| s /M string ~ M GU T/(2e 3 ) ~ 7 x 10 14 GeV, which is just right for 
generating acceptable light neutrino masses via the see-saw mechanism. Although at the present 
operator order the Majorana mass terms vanish with the non-vanishing vevs discussed earlier, 
non-trivial operators may exist at higher order. 

V. CONCLUSION 

In this paper we construct three-family PS models in the Zg abelian symmetric orbifold. Our 
models are mainly motivated by recent discussions on orbifold GUTs. We are able to realize some 
features of the orbifold GUTs in the string compactification limit where the compactified space is 
effectively 5d. The breaking of the N = 2 to N = 1 supersymmetry in 4d and the Eg (or SOio) 
gauge symmetry to that of PS are realized similarly as in orbifold GUTs. We find three family 
chiral matter fields, two of them can be regarded as "brane" states and one as a "bulk" state, in the 
terminology of orbifold GUTs. These models extend three-family orbifold string model building to 
non-prime-order orbifolds, and we find some new features in the matter spectra when compared 
to that of the prime-order orbifolds. Matter fields arise not only from the untwisted but also from 
twisted sectors, and typically there is a horizontal 2 + 1 splitting in the family space. This splitting 
may have the potential to better facilitate the description of fermion masses and mixings. 

We find one of our string models, with an E6 gauge group in 5d, is particularly interesting. It 
has the following properties: 

1. Renormalizable Yukawa couplings exist only for the third family. Moreover the model pre- 
dicts a unification relation among the third family Yukawa couplings and GUT gauge cou- 
pling at the 5d compactification scale. 

2. The renormalizable and non-renormalizable couplings can affect a spontaneous breaking 
of the PS symmetry to that of the SM with fields in the (6, 1, 1) and (4, 1, 2) + (4, 1, 2) 
representations. Moreover, after the symmetry breaking, many unwanted states can develop 
large mass. 
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3. There is a non-abelian D4 family symmetry for the first two families, which makes the 
model a good playground for studying fermion mass hierarchy and the flavor problem in 
supersymmetry. 

We regard these phenomenological features as merits of our models. 

On the other hand, there are two main problems for the models (of course, they are not unique 
to our models). 

• Exotics: There is a pair of SM vector-like exotic particles at the low energy scale. To the 
order of our analysis, we have not found any Yukawa coupling that can give them mass. 

• Proton stability: We have Higgs fields transforming in the same (or conjugate) PS represen- 
tation as the SM right-handed matter. In general, baryon/lepton- number violating effective 
operators are induced after PS symmetry breaking. We have not yet found a symmetry that 
can distinguish the PS breaking Higgses from matter and effectively eliminate the dangerous 
operators to all orders. 

These problems have only been examined at a rather primitive and qualitative level. Specifically 
we have looked for non-trivial Yukawa couplings allowed by string selection rules to certain orders. 
We have shown that the baryon/lepton- number violating effective operators can be avoided if 
one prudently chooses appropriate values for the singlet vevs. Unfortunately, we are not able to 
extend this argument to all orders in string perturbation. The problem apparently results from 
the presence of only the trivial fixed point in the T\ twisted sector on the G2 torus. This has the 
consequence of nullifying any space group selection rules for the G2 torus whenever any T\ twisted 
sector operator is present. 

In any case, these problems may point to the direction for refining our models. For example, it 
would be desirable to search for models with fewer color triplets and, if possible, no exotics. The 
analysis presented here is clearly just the beginning. It would certainly be useful to expand the 
search to other Z2 x Z^r orbifolds with one (or two) Wilson lines in the SO4 direction in order to 
find more effective 5 or 6d orbifold GUTs. (The Z2 x Z2 model might be particularly interesting 
and it has been studied recently in ref. [2^- Presumably the model is simpler than ours because 
it has only three twisted sectors and the number of modular invariant gauge embeddings is more 
limited. However, realistic three-family model have yet to be constructed.) In brief, we believe our 
analysis has opened up a promising new direction for string model building. 
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APPENDIX A: RECIPES FOR CONSTRUCTING NON-PRIME-ORDER ORBIFOLD 
MODELS IN THE PRESENCE OF DISCRETE WILSON LINES 



In this appendix we review the construction of non-prime-order orbifold models with discrete 
Wilson lines. 

Our starting point is the lOd heterotic string theory, which consists of a 26d left-moving bosonic 
string and a lOd right-moving superstring. Modular invariance requires the momenta of the internal 
left-moving bosonic degrees of freedom (16 of them) lie in a 16d Euclidean even self-dual lattice, 
we choose to be the Es x Es root lattice. 12 

To make a connection to the 4d world, we must compactify 6 spatial dimensions of the 10 
remaining space-time dimensions. There are many ways (see, e.g., ref. [j| for a collection of 
early works) to achieve a 4d N = 1 supersymmetric spectrum, among them the most studied 
in the literature is the orbifold construction 



BE 



10 



11 



. We will use the simplest abelian 



symmetric orbifold construction, where we differentiate the space-time and internal degrees of 
freedom, and realize the orbifold twists and Wilson lines by shifts in the Es x Es lattice. This type 
of construction admits a clear space-time interpretation. 

The full definition of an orbifold model requires the specification of a six-torus T 6 , corresponding 
to the compactified spatial dimensions, a point group, (such as the cyclic groups Zjy or Zjy x Zjvf 
with N, M = 3,4,6,7,8,12 Q), corresponding to the automorphism of the T 6 -lattice, and an 
embedding of the space group, (which consists of the point group and the lattice translations), 



12 For an orthonormal basis, the Eg root lattice consists of following vectors, (m, ri2, ■ ■ ■ , n&) and (ni + \,n% + 



i • • ■ , ng + 4), where m, ri2, ■ ■ ■ ng are integers and X!i=i n i ~ m °d 2. 
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in the Eg x Eg lattice. Normally one denotes the generator of the discrete group Zjy by a twist 
vector v = (v±, V2,vs), acting on the three complex planes by 6 : Z{ — ► e 27T1Vi Zi (i = 1,2,3). 
To ensure that one space-time supersymmetry survives in 4d, the twist vector needs to satisfy 
±v\ ± V2 ± V3 = mod 2, 13 and none of the v^s vanishes. In the abelian orbifold construction, 
embeddings of the space group in the Eg x Eg are realized by shifts of the corresponding lattice, 
P — > P + k\ + ZW, where k, I are integers, P are vectors of the Eg x Eg root lattice, V the gauge 
twists, realizing the point group, and W the discrete Wilson lines, realizing the lattice translations. 
The cyclic group multiplication rules require A^V and Nyy W to be in the Eg xEg lattice. (In general 
the degrees of the Wilson lines, Nw, divide the degree of the orbifold twist, N.) 

String states closed on T 6 , i.e. those satisfying the condition Zj(r, a + ir) = Z{(t, a) modulo lat- 
tice translations, give rise to the untwisted-sector states. Besides the N = 1 supergravity multiplet 
and modulus fields that parameterize deformations of the background fields, the untwisted-sector 
states also give rise to gauge and matter fields. Embeddings of the point group in Eg x Eg break 
the gauge symmetry down to its commutator subgroups, i.e., the surviving non-zero roots satisfy 

P 2 = 2, P-VeZ. (Al) 

Note that we cannot lower the rank of the surviving groups in abelian orbifold models, since by 
construction the gauge twists and Wilson lines commute with the Eg x Eg Cartan subalgebra. 

The conditions for the untwisted-sector matter states are similar. It is convenient to bosonize 
the right-moving fermionic degrees of freedom and denote their SOg momenta in the light-cone 
gauge by r. 14 The right-moving NS sectors (&!_]/ 2 f° r i = 1:2,3) of the untwisted matter states 
have SOg weights ri = (0,1,0,0), T2 = (0,0,1,0) and r3 = (0,0,0,1) and they pick up phases 
e -27nivv unc j er the orbifold twists, so the corresponding roots must satisfy 

P 2 = 2, P-V-r r veZ (i = 1,2,3). (A2) 

This gives three untwisted-matter sectors Ui, one for each complex plane. (Eg. lAll can also be 
written in the same form with ro = (1, 0, 0, 0), where the non-zero entry lies in the uncompactified 
direction) . 

The gauge groups and untwisted-sector matter spectra are modified when discrete Wilson lines 
are turned on. In the presence of general background fields Gij, Bij and Wj, the canonical momenta 

13 The signs are arbitrary. We will use the convention that all signs are positive. 

14 That is, r's are in the SOs weight lattice, the integral and half-integral weights correspond to the Neveu-Schwarz 
(NS) and Ramond (R) sector fermions. For R-sector weights, the first components, ro, indicate the helicities of 
space-time fermions. In this notation, the first component of the twist vector, vo, is zero, r are commonly referred 
to as the H momenta. 



32 



conjugate to the compactified coordinates and the gauge coordinates are IT* = p\ + p R + \ W l • 
(Il + p L ) + B l i(p L j —Pjij) and II = p L — ^Wi(p L —p R ). Since v\ — Vr = 2n l , = m l and II = P, 
where the integers m ! ,n ! are the momentum (or KK modes) and winding quantum numbers, we 
(the string unit is a' = 1/2) 



have 
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PL = —+(G^-B^--W^W^n j --P-W\ 

p R = pi - 2n\ p L = P + n'Wi . (A3) 

In the models studied in this paper, we take Bij = and Gij = ^RiRjei ■ ej is defined by the 
geometry of the internal space T 6 with basis vectors given by ej and dimensions by Ri. In this 
case, these equations, along with the masslessness conditions 

-va\ = N R + \G ijP y R + ir 2 - i = (A4) 

\m 2 L = N L + ±G ijP y L + ip| - 1 = (A5) 

where Nl, Nr are integral oscillator mode numbers and the last two terms in eq. IA4I are the 
contribution of the bosonized NSR fermions, require the winding number for massless states to be 
zero, and the last term of p l L an integer (for both gauge and matter fields), i.e., 

P • W € Z . (A6) 

String states closed on themselves under the identification of a non-trivial element of the 
point group (modulo translations by lattice vectors) give rise to the twisted-sector states. For 
the fc th twisted-sector (for which the complex compactified coordinates satisfy Z{(r,a + ir) = 
e 2mkVi Zi(T,a)), the Es x Es and SOs momenta are shifted according to P — > P + k~K n( and 
r — > r + Arv, where X nf = V + nJ-Wj with n\ being fixed-point dependent winding numbers. 15 The 



massless states satisfy the following equations 



1 6 



Y,NM k) + - 2 (r + kvf + a^=0, (A7) 

i=l 

\m\ = f: NfrW + \(P + kX nt f + af = , (A8) 
i=i 

where Nf' and are intergral numbers of the right- and left-moving (bosonic) oscillators, a R , 
a £ are the normal ordering constants, 



4^ = -« + l^2\ kv i\ f 1 - \ kv i 



i=l 



15 Discrete Wilson lines break the degeneracies of the fixed poin ts. The integers n\ should be chosen appropriately, 
depending on the direction and degree of the Wilson line |25l . li^ . . 
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1 3 

{k) - -l + ^l^il (l-Nil) , (A9) 

i=l 



L 2 



with kv i = kv-i mod 1, and u;- = kvi if kv. t > and 1 — kvi if kv-i < are oscillator energies. Thus, 
a twisted-sector matter state can be labelled by the twisted sector, k, the fixed point it is localized 
on, f , the number of windings associated with the fixed point, nj. , and the number of right- and 
left-moving oscillators, N? 1 and Nf. 

Note, however, that not all gauge twists and discrete Wilson lines are physically allowed. To 



ensure modular invariance of the one-loop partition 
the right- and left-movers, one needs to require P, |s 



function (or the level-matching condition) for 



10 



24. 



N(X. 2 nf - v 2 ) = 0mod2. (A10) 

In addition, in non-prime-order orbifold models, the degrees of Wilson lines are in general the di- 
visors of that of the orbifold twist; they need to satisfy more stringent modular- invariance require- 
ments. For example, in the model where v@ = g(l,2, — 3) there are at most three admissible 
Wilson lines, two of degree-2 (W^, i = 1,2) and one of degree-3 (W3) 



43, |4^. 16 Modular- 



invariance conditions involving these Wilson lines are the following more restrictive set (where 
* = 1,2), 

{2(W^) 2 , 3(W 3 ) 2 , 4W^ 1} • W^ 2) , 12W^ } • W 3 } = 0mod2. (All) 

Further complications arise for non-prime-order orbifold models. Unlike the prime-order orb- 
ifolds such as the Z3 and Z7 models, fixed points of the higher- twisted sectors are not always 
invariant under the defining orbifold twist, 9. However, one can find linear combinations of the 
states corresponding to these fixed points such that they have definite eigenvalues under the 9 
rotation. As it has been shown in ref. [313; for any fixed point represented by a space-group 
element (9 k ,l) (i.e., a fixed point f satisfying f = 9 k f + 1, where 1 is a vector of the T 6 lattice and 
called the equivalent shift vector), if it can be written as a power of a prime element (0 m ,Y) (a 



16 Actually, the number of admissible Wilson lines also depends on the compactified lattice. If we restrict it to a Lie 
algebra root lattice, then there are four possibilities, (A) SU 6 ffi SU 2 , (B) SU 3 © S0 8 , (C) SU3 © SO7 © SU 2 , and 
(D) G2 ffl SU3 ffi SO4, whose Coxeter elements realize the orbifolding. (Note that the lattices C and D are the 
sublattices of B, where the quotients are isomorphic to Z2 and Z2 x Z2. The Z2 x Z2 is the center of SOs- There 
are two additional lattices, SU4' 2 ' © SU3 ffi SU2 and SU3' 2 ' © SU3 ffi SO4, which also involve outer automorphism of 
the root lattice and give identical results to C and D.) The A (B/C) lattice has at most one degree-2 (one degree-2 
and one degree-3) Wilson line(s). The D lattice corresponds to the case in the main text and is the most intuitive 
one. For our three-family models, we choose to turn on two Wilson lines, one degree 2 and one degree 3. The low 
energy phenomenology for the last three lattices are different, since both the quantum numbers (which label the 
states) and the selection rules (which determine allowed couplings) depend on the lattice choice. In particular, the 
D4 family symmetry does not subsist for the B and C lattices. 
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prime element is the one that cannot be written as a power of any other element), then the linear 
combinations, 

m— 1 

\k,l) = ^j- e \9 k ,9h), (A12) 
e=o 

have eigenvalues 7 = e 27rm / m and n = 0, 1, • • • , m — 1 under the # rotation. We can therefore use this 
notation to appropriately label twisted-sector states arising from different fixed points according 
to their ^-eigenvalues. 

Finally, in the Zq orbifold that we are most interested in (and other orbifolds containing sub- 
orbifolds with fixed tori), special care must also be taken in the presence of discrete Wilson lines. 
This orbifold contains Z2 and Z3 sub-orbifolds with twists {1,# 3 } and {1,# 2 ,# 4 } respectively. The 
second and third complex planes are invariant tori of the 9 3 and 6 2 /6 4: twists, i.e., they are unrotated 
under the respective orbifold twists. Consequently these directions have mode expansions of a 
toroidal coordinate in the T3 and sectors, as in eq. IA31 with gauge momenta replaced by the 
shifted momenta, P + 3X„ f in the T3 and P + 2X nf , P + 4X„ f in the T 2 , T4 sectors. There are also 
Wilson- line dependent terms. The mode expansions of the second direction in the T3 sector and the 
third direction in the T 2 ^ sectors contain the 2W3 and 3W 2 ^ (i = 1,2) Wilson lines respectively. 
(The numeric factors are due to the fact that the only non-trivial action for the W3 and W 2 ^ 
Wilson lines are the Z2 and Z3 sub-orbifold twists.) Just like the untwisted-sector states, the 
winding numbers for massless states along these unrotated directions in respective twisted sectors 
must be set to zero, and the shifted momenta must satisfy the (necessary) projection conditions 
2(P + 3X nf ) • W 3 e Z for the T 3 and 3(P + 2X„ f ) • G Z, 3(P + 4X„ f ) • G Z for the T 2 , T A 
sector states. Using the modular invariance conditions, eqs. lAlOl and lATTl these conditions reduce 
to 

P • W3 G Z , for T 3 sector ; P • wf ) G Z , for T 2A sector . (A13) 

These projections are crucial for rendering an anomaly-free mass spectrum and have not been 
properly addressed in the literature. 

Multiplicities of the twisted-sector states are computed from the generalized Gliozzi-Scherk- 
Olive (GSO) projector [lfll . which can be derived from the one-loop partition function by 
power-series expansions. In our notation, for a twisted-sector state labelled by {fc,7,nf} and 
appropriate oscillator numbers, the projector is [4^] 



N-l 

^,7,n f ) = -^[A(fc, 7 ,n f )f , (A14) 



1=0 
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where 



A(k, 7, nf) = 07exp < in 



(2P + fcX nf )-X nf -(2r + A;v)-v 



(A15) 



In this expression, cj> = exp[2-7ri Ylt=ii^t ~ Nf~)vi] is the oscillator phase, with V{ = sga(kvi)vi, 
= — sgn(kvi)vi if kvi 7^ 0, and Vi, Vi+3 = if kvi = 0, for i = 1,2,3. Only states with 
A(fc, 7, nf) = 1 will survive the projection. (Equivalently, the generalized GSO projector can be 
defined as in ref. ^^], in a slightly different fashion. We find it is more convenient to implement 
eq. IA14I in the non-prime-order orbifolds.) 

The GSO projector in eas. lA14l and lA15l needs to be modified slightly for states in the untwisted 
sectors and those twisted sectors with fixed tori. In the untwisted case where k = 0, 7 = = 1, 
nf = 0, the required modification is simply an additional projection with respect to the Wilson 
line, 



A' 



N w -1 

Nw £ exp 



mf =0 

P W € 



ivr^(2P • X mf - 2r • v) 



(A16) 



the same projection conditions as in eqs. IA11 IA2I 



Then A = 1 gives P • V — r • v G 
andlA"6l 

The second case is more complicated. Assume that the Tf. sector has a fixed torus, and the 
Wilson line of degree N\y also lies in this torus. Then we must have N\y divides k, and the required 
modification to the GSO projector is 

^ N w -1 



A ( 



N w 



exp < m£ 



mf =0 



(2P + kX mt ) ■ X mf -(2r + fcv) • v) 



(A17) 



In this expression, it must be understood that each shifted momentum, P + /cX mf , satisfies its own 
masslessness condition, eq. IA8I Therefore the momenta P in different terms differ by a multiple 
of fcW, which is a vector in the Es x Es root lattice; these states are isomorphic to each other. 
Requiring that A = 1 in ea, IA17l then implies an additional projection P-W G Z for the T/% twisted 
sector states. It can also be seen that with this condition eq. IA17I reduces to eqs. IA14I and IA15I 
with X nf = V. These discussions can be easily generalized to the cases with several Wilson lines. 



In the Zg model, they give the same projection conditions as in eq. IA13I 

While it is quite tedious to carry out all these procedures in practice, they can be easily com- 
puterized, as follows. First we check the modular invariance of the gauge twists and Wilson lines, 
and find all the Es x Es roots that satisfy eqs. lAll and IA6I From these we find the set of linearly 
independent positive roots (i.e., the simple roots) and determine the unbroken gauge groups. Next 
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we find the roots corresponding to the untwisted- and twisted-sector matter states, from eqs. IA21 
IA6I and eqs. I A 71 IA8I respectively. For the twisted-sector matter, we only keep states of negative 
helicity (i.e. left-handed), and compute their multiplicities using the generalized GSO projector in 
eq. IA141 We also need to take special care of the twisted-sector states in the presence of Wilson 
lines, as noted in the paragraph after eq. IA12I We then decompose the roots into highest weight 
representations, and determine their Dynkin indices and dimensions using the standard group the- 
oretical method. Finally, the remaining abelian factors are found by searching for orthogonal roots, 
and the charges for matter states are determined accordingly by projecting the matter states onto 
these abelian roots. 



APPENDIX B: THREE-FAMILY PATI-SALAM MODELS 



1. The models 



In this appendix, we apply the procedure outlined in appendix to construct three- family PS 
models in the 7<q orbifold with V6 = 2, —3). The Zg is equivalent to the %i x Z3 orbifold, where 
the two twist vectors are V2 = 3vg = 0, —1) and V3 = 2v@ = |(1, —1,0). 

There are in total 61 inequivalent modular invariant choices for the gauge twists in the 7jq orb- 
ifold model [2J]. To narrow down the possibilities, we demand the models we start with (before 
imposing any Wilson line) contain an SO10 gauge group and some matter fields in 16/16 repre- 
sentations in the first or third twisted sectors. Although this step makes our results less generic, 
it greatly reduces the large number of possible models to a manageable subset. We choose the 
following two gauge twists, 

• model A: 

V 6 = - (22200000) (11000000) , (Bl) 
6 

• model B: 

V 6 = - (41100000) (22000000) , (B2) 
6 

which break the Eg x Eg gauge symmetry down to SO10 x SU3 x E-j' and SO10 x SU2 x E7' respectively. 
The model A (B) contains four (two) 16 and one (three) 16 in the untwisted sectors, and eighteen 
(fourteen) 16 and three (seven) 16 in the twisted sectors; in total there are eighteen (six) SO10 
families 12411 . 
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To further break the gauge symmetries and reduce the number of families, we impose discrete 
Wilson lines. As previously mentioned, there are at most one degree-3 Wilson line in the second 
complex plane, and two degree-2 Wilson lines in the third. We choose to add two of them, one of 
degree- 2 and one of degree-3, as follows, 

• model Al: 

W 2 = \ (10000111) (00000000) , W 3 = \ (1 - 1000000) (00200000) , (B3) 

• model A2: 

W 2 = - (10000111) (00000000) , W 3 = - (21 - 100000) (02110000) , (B4) 

• model B: 

W 2 = ]- (00011000) (10100000) , W 3 = - (01 - 100000) (00200000) . (B5) 

One can easily verify our choices satisfy the modular-invariance requirements, eqs. lAlOl and IA11I 
The W3 Wilson line of models Al /B and A2 breaks the gauge group in the observable sector to 
SO10 and SO10 x SU2, respectively, and the W 2 breaks them further down to the PS gauge group, 
in all three models. 

The remaining unbroken gauge groups are SU4C X SU 2 l X SU 2 r x SOio' x SU 2 ' x (U1) 5 for 
models A1/A2 and SU4C x SU 2 l x SU 2 r x SOio' x (U1) 6 for model B. The untwisted- and twisted- 
sector matter provide the following irreducible representations of the PS gauge group (modulo 
some singlets), 

• model Al: 

U X '. (4,2,1), U 2 : (1,2,2), U 3 : (4, 1, 2) + (4, 1, 2) , 

T\ : 2(4, 2, 1) + 2(4, 1, 2) + 4(4, 1, 1) + 4(4, 1, 1) + 8(1, 2, 1) + 8(1, 1, 2) + 2(1, 1, 2; 1, 2) , 
T 2 : 2(4, 1,2) + (6, 1,1), T 3 : 6(6, 1, 1) + 6(1,2,2) , T 4 : (4, 1, 2) + 2(6, 1, 1) , (B6) 

• model A2: 

Ui : (4, 2, 1) , 

T\ : 2(4, 2, 1) + 2(4, 1, 2) + 4(4, 1, 1) + 4(4, 1, 1) + 8(1, 2, 1) + 6(1, 1, 2) 

+ 2(1, 2,1; 1,2) + 2(1, 1,2; 1,2), 

T 2 : 2(4, 1,2) + (6, 1,1) + (1,2, 2), T 3 : 2(4, 1, 1) + 2(4, 1, 1) + 6(1, 1, 2) , 

T 4 : (4, 1,2) +2(6, 1,1) +2(1, 2, 2), (B7) 
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• model B: 

U 1 : (4,2,1), U 2 : (6,1,1), U 3 : (4, 2, 1) + (4, 2, 1) , 

7i : 2(6, 1, 1) + 2(1, 2, 2) + 4(4, 1, 1) + 4(4, 1, 1) + 10(1, 2, 1) + 8(1, 1, 2) , 

T 2 : 2(4, 1, 2) + (1, 2, 2) , T 3 : 4(4, 2, 1) + 4(4, 1, 2) + 2(4, 2, 1) + 2(4, 1, 2) + 6(1, 2, 1) , 

T 4 : (4, 1,2) +2(1, 2, 2). (B8) 

In general, the T 6 _fc-sector is the CPT conjugate of the T^-sector, therefore T5 does not give rise to 
additional states. However we need to keep those states from the T4-sector whose CPT conjugates 
in the T2-sector have positive helicities. 

We see there are indeed three chiral families ((4, 2, 1) + (4, 1, 2) under the PS group), two from 
the T\ (T3) twisted sector and one from the untwisted and the T 2j 4 twisted sectors in model A1/A2 
(B), modulo some (4, 1, 2) + (4, 1, 2) (and (4, 2, 1) + (4, 2, 1) in model B) vector-like pairs. Each 
(4,2,1) + (4,1,2) family encompasses one complete family of the SM quarks and leptons (plus 
the right-handed neutrino), since they decompose under the SM gauge group, SU3C x SU2L x Uiy, 
as follows: (4, 2, 1) = (3, 2) 1/6 + (1, 2)_ 1/2 , (4, 1, 2) = (3, l) 1/3 + (3, l)_ 2/3 + (1, l)j + (1, 1)„. 
The complete matter spectra for these models is listed in appendix lB 21 where the notation for the 
twisted-sector states is also explained. Three-family PS models have been previously constructed 
in heterotic string models using the free-fermionic method but to our knowledge, has not been 
realized with abelian orbifolds. 

Like all the three-family orbifold models in the literature , LU| , our models suffer from 



the embarrassment-of-riches syndrome, i.e., they contain many candidates for SM exotic particles, 

(4.1.1) , (4,1,1) and (1,2,1), (1,1,2). (These states decompose under the SM gauge group as 
follows: (4,1,1) = (3,l) 1/6 + (l,l)_ 1/2 , (4,1,1) = (3, l)_ 1/6 + (1, l) 1/2 , (1,2,1) = (1,2) and 

(1.1.2) = (1,1)1/2 + (1)1)— 1/2) an d have exotic charge assignments.) Nevertheless, these exotic 
states are necessary for the quantum consistency of the models. Models A1/A2 (B) contain five 
(six) additional abelian symmetries (their charges for the matter fields are listed in appendix lB 2j) . 
one of them is anomalous. This Uia anomaly satisfies the following universal condition, 



TiT(R)Qa = TtQIQa = ^TrQi = ^TrQ A 



■^g, Model Al, 



Model A2 , (B9) 



§ , Model B , 



~ 1/2 

where Qi t A = Qi,A/k i A with ki,k A are the normalization factors ("levels") of the corresponding 
abelian groups, and 2T(R) is the index of the representation R under a specific non-abelian factor. 
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This condition guarantees the anomalies are cancelled by the generalized Green-Schwarz mechanism 
[4^ in 4d 43]. In fact, considerations of these anomalies provide useful consistency checks for our 
computer program. One can also check that all the 4d chiral anomalies for the non-abelian gauge 
groups and global anomalies for the SU2 groups vanish. 

We note that our models differ from the three-family SM-like models constructed from Z3 
orbifolds in several respects. For example, in the simplest model of the first paper in ref. 
the left-handed quarks arise solely from the untwisted sectors, one family each from the three 
untwisted sectors. The right-handed quarks and left-handed leptons arise from the twisted sectors, 
and the multiplicities of three come from the three equivalent fixed points in one of the complex 
planes. As such, the three families are completely degenerate in the horizontal family space. Our 
models exhibit quite different spectral patterns. Two complete families come from the T1/T3 
twisted sectors, and one family comes from the combination of the untwisted and T2/T4 twisted 
sectors. This type of pattern breaks the degeneracies among the three families, and may have 
better prospects for explaining the observed fermion mass hierarchy and CP phases. It will also 
be interesting to generalize our models to other non-prime-order and Z^r x %m orbifolds to see 
whether this feature persists. 

The above three-family PS models have some advantages over string-based SM-like and GUT 
models which have been the foci of model buildings in the past. The SM-like models usually contain 
many extra abelian gauge symmetries, and in general it is difficult to reproduce the standard 
hypercharge normalization as in GUT theories. (For an extensive search of Uiy in the Z3 orbifold, 
see ref. [5(1].) This may upset the successful prediction of the weak- mixing angle in the context of 
supersymmetric unification theories (see, however, ref. for alternative tri- unification models). 
In the PS model, on the other hand, the hypercharge normalization is standard, by construction. 
The hypercharge Uiy is a diagonal subgroup of the Uit 3 C SU2R and the U 1 ( B _ L ) C SU4C, i.e. 
Qy = T3R + \(B — L), hence /cy = l + |x| = |. Moreover, all our models contain SO10 spinor 



representations, which may be a welcoming ingredient for better low energy phenomenology [52 1. 
(In contrast, the spinor representations are absent in SM-like string models and models based 
on the Spin(32)/Z2 heterotic string and superstrings.) Finally, when compared with string GUT 
models based on SU5 or other larger gauge symmetries jl^T ]. the PS models generically do not need 
large Higgs representations (such as those in the adjoint representation) to break the GUT gauge 
symmetries. These large Higgs representations require one to realize the current algebra of the 
gauge symmetries at higher Kac-Moody levels jj^J and render the models more complicated. 
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2. Complete matter spectra 



We list the complete matter content of the three-family PS models, including all the singlet 
fields, in tables ITU ITTTI and HVI 

Let us explain our notation for the twisted-sector states somewhat. For the sake of presen- 
tation, we take the six compactified dimensions to be a factorizable Lie algebra root lattice, 
G2 © SU3 © SO4 (see fig. EJ), as in sect. 1111 Al 17 The notation, although pertaining to the par- 
ticular case, can be generalized to any other lattice. 

The twisted sectors contain various matter states, arising from different fixed points. It is useful 
to introduce the following notation for the twisted-sector states (see appendix also ref. 4a] h 



T Hl) a ("-3)("-2n' 2 )N, (BIO) 

where the first sub-index represents the k th twisted sector, the next three quantum numbers, (7)^, 
(713) and (77,2, n 2 ), specify fixed points on the first, second and third complex planes, respectively. 
The last sub-index ./V = (NiN 2 N^) denotes an array of three integral oscillator numbers for the left 
movers. For example, (100) and (100) represent creating one oscillator in the Z\ and Z\ directions. 
We only need to keep the left-moving oscillator numbers since all the massless states in our models 
have zero numbers of right-moving oscillators. 

The quantum numbers (7)0, (77,3) and (712, n 2 ) label conjugacy classes of the fixed points. The 
second complex plane, i.e., the SU3 lattice, is a Z3 orbifold, and has three fixed points, f = 0, 
q-(2e3 + 64) and g(e3 + 2e4), where e3, are the basis of the SU3 lattice. The conjugacy classes 
are labelled by equivalent shift vectors ae^ + be^, where + 6 = 713 = 0, 1,2 mod 3 is the "winding 
number." We can introduce at most one degree-3 Wilson line, and specify these fixed points by 
723 = 0,1,2 for k 7^ 3, and 77,3 = for k = 3 (since Z 2 is the fixed plane of the 6 3 twist). 
Similarly, the third complex plane, i.e., the SO4 lattice, is a 7L 2 orbifold, and has four fixed points, 
f = \{n 2 e^ + n' 2 eo), where n 2 ,n 2 = 0, 1 and es, e$ are the basis vectors of the SO4 lattice. The 
equivalent shift vectors for these fixed points are —n 2 es — n 2 eg. We can introduce at most two 



17 In fact, the matter spectra also depend on our choice of the compactifed lattice. The Lie algebra root lattice sits at 
special points of the six-torus modulus space with enhanced gauge symmetries. The spectra (if we restrict ourselves 
to the root lattice) should be regarded as truncated spectra of the complete model, where there are additional 
states neutral under the observable gauge groups. The use of non-simply-laced algebra G2 may also be worrisome. 
But it may be constructed from a simply-laced algebra such as the SOs by appropriate outer-automorphism, at the 
expense of reducing the rank of the hidden-sector groups. The G2 lattice can also be replaced by the SU3' 2 ' lattice 
without changing any of our discussions on low energy phenomenology. The SU3' 2 ' Coxeter element is generated 

by the Weyl reflections with respect to ei, Si, and with respect to e 2 — ei, S2<->i, Cq [2] = siS2<-»i 



-1 1 



where ei and e2 are the simple roots of the SU3. 
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TABLE II: Matter spectrum of model Al. The levels of the Ui groups are 24, 24, 72, 36 and 48. 
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TABLE III: Matter spectrum of model A2. The levels of the Ui groups are 16, 72, 32, 80 and 36. 
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TABLE IV: Matter spectrum of model B. The levels of the Ui groups are 24, 88, 40, 216, 88 and 576. 
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independent degree-2 Wilson lines, and specify these fixed points by a pair of winding numbers, 
n2,n' 2 = 0, 1 for k = 1,3,5, and n-i = n 2 = for k = 2,4 (since Z3 is the fixed plane of the 8 2, 
twists). In the present models, we introduce only one non- vanishing Z2 Wilson line, so the degree-2 
degeneracy associated with n 2 is not resolved. 

For the first complex plane (i.e., the G2 lattice), there is only one fixed point in the T\ twisted 
sector, f = 0. The corresponding state will be denoted by (7) = (1). In the T 2 and T4 twisted 
sectors, there are three classes of fixed points, f a = |aei with a = 0,1,2. The last two, fi^, 
transform into each other under the orbifold twist 8, while the first, fo, is invariant. For states 
associated with the non-invariant fixed points, we can take linear combinations such that they 
have definite ^-eigenvalues, 7 = ±1. We label the two eigen-states with 7 = 1 by | fo) = (l)i, 
-^(|fi) + I F2) ) = (1)2 and the one with 7 = — 1 by ^(|fi) — | fa) ) = ( — !)• Similarly, there are four 
classes of fixed points in the T 3 twisted sector, f a & = |(aei + be2), with a, b = 0, 1. Three of them, 
foi, fio and fn, transform into each other under the orbifold twist 8, while the other, foo, is invariant. 
After taking linear combinations, we have two states with 7 = 1, | foo) and ^=(|fio) + |foi) + |fn)), 
denoted by (1)1,2, and two states 4g(|fro) + 7|foi) + 7 2 |fn)) with eigenvalues 7 = 0;, u> 2 , where 
to = e l27r / 3 , denoted by (u) and (u> 2 ). 

Multiplicities of the twisted-sector states are computed from the generalized GSO projector of 
eq. IA141 They equal the number of degenerate fixed points with quantum numbers k, 7, a, 713, 
ri2, n 2 and N. With two Wilson lines, one degree-2 and one degree-3, the Ti-sector states have 
multiplicities two, due to the unresolved degree-2 degeneracies in n 2 . The T2 and T4-sector states 
have multiplicities two and one, for states with eigenvalues 7 = 1 and 7 = — 1, since they have two 
and one associated fixed points, (l)i,2 and (—1). Finally, the T3-sector states have multiplicities 
four and two, for states with eigenvalues 7 = 1 and 7 = lu,uj 2 , since the associated fixed points are 
(l)i,2 and (w), ipj 2 ) and there are additional degeneracies in n 2 = 0, 1. 

APPENDIX C: YUKAWA COUPLINGS IN THE Z 6 MODELS 
1. Selection rules 

In this subsection, we give the complete set of string selection rules in the Zq orbifold model. 
Although these rules are well known, they are not always stated correctly in the previous literature. 
In heterotic string models, physical states corresponding to space-time bosons can be written 
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in terms of vertex operators in the (— l)-ghost picture, they are 

V W = e-^n ( aZ0^ ) (^)^ ) e- 2irW - H e 2ip(f) - x af, (CI) 

i=l 

where cj) and (dZi) N i \ (dZ{) Ni denote the superconformal ghost and and N^p left-moving 
oscillators in the Zi and Z% directions. P^ and are the Es x Es and SOs shifted momenta (X 
and H are bosonic coordinates parameterizing the respective maximal tori); the latter is commonly 
known as the H-momentum. More explicitly, in the model the internal SC-6 parts of the H- 
momenta for states in different sectors are: 

Ui: (1,0,0), U 2 : (0,1,0), U 3 : (0,0,1), 

Tf. (l,2,3)/6, T 2 : (2,4,0) /6, T 3 : (3, 0, 3)/6, T 4 : (4, 2, 0)/6. (C2) 

Finally, o j denotes the twisted field, creating twisted-sector vacua out of the untwisted ones; it 
is thus trivial for the untwisted sector states. 

Equivalently, due to the world-sheet superconformal symmetry for the right-moving superstring, 
bosonic state vertex operators can be written in other ghost pictures by acting with a picture 
changing operator Q, D e^^LiWi^i + $$Zi) = e^ £Li( e2ir "' H ^ + e" 2 <- H dZ;), where 
r\ = (1,0,0), r 2 = (0,1,0), = (0,0,1). Therefore in the 0-ghost picture the H momenta 
are reduced by r^'s and additional factors of dZi and dZi are introduced. However the unique 
combination (defined up to a modding with respect to the degree of orbifolding) 

flf =r< W -JVf (C3) 

remains invariant under picture-changing operations. 

For Yukawa couplings in string theory, one can use the standard conformal field theory technique 
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551 ] to compute the corresponding ?i-point correlation function, 



{V^%V%V^---vt\ (C4) 

where Vila's are the vertex operators for space-time fermions in the (— l/2)-ghost picture and all 
but one of the bosonic vertex operators have been brought into the 0-ghost picture. Note that the 
total ghost charge must cancel with a background ghost charge of 2. 

For our purposes, however, we do not need the exact form of the above correlation function; 18 
string selection rules suffice to tell us whether certain types of Yukawa couplings are non-trivial, 
i.e., not identically zero. These rules are provided by different parts of the vertex operator. 



The functional dependence of Yukawa couplings on moduli can be determined by a field theoretical method 
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The gauge part requires conservation of the Es x Es momenta, = 0, which is nothing but 

the gauge invariance condition for matter couplings. A similar rule also holds for the H-momenta, 
but we have to be more careful because these momenta are not invariant under the picture-changing 
operations. Instead of the H-momenta, the rules must be imposed on the well-defined R-charges 
of eq. IC31 they are 

R? = 1 mod 6, R 2 ] = 1 mod 3 ' Yl R s ] = 1 mod 2 - ( C5 ) 

lie. 

We note that they can be understood as conservations of discrete global R-charges in the field 
theory language 60] . The nature of R-symmetries can be seen from the spacetime supersymmetry 
charge Qsusy = / ^e _< ^/ 2 5'e 2ir susY' H where S is the spin field for uncompactified dimensions, 
r SUSY = ^(1? 1 ; 1) 3 and the integration is over the world-sheet coordinates. Under the orbifold 
twist Qsusy — ► e^QsusY and the bosonic vertex operator — > e~ 2l7rR ; Vi V^. Eq. IC5I is 
equivalent to the statement that the superpotential is invariant under this twist. 

The part associated with the twisted fields leads to the so-called space group selection rule. 
Global monodromy requires that the product of space group elements associated with the fixed 
points contain the unit (1,0), that is, 

^A^ = 0mod6, (C6) 

i 

e i 

where the indices t sum over the twisted-sector states (in the &^-th sector), denote the fixed 
points under these twists, and are arbitrary lattice vectors of the compactified space. Eq. IC6I 
is the point group selection rule for the Zg models. 

Eq. IC7I implicitly depends on the compactified lattice. For our choice of the root lattice 
G2 © SU3 © SO4, using the notation of appendix IB 21 we find it is equivalent to 4|| 

Y = mod 3 , ^(n^, nf ] ) = (0, 0) mod 2 , (C8) 



The low energy effective action of 4d strings enjoys a target space duality symmetry PSL(2,Z) [5^ | (or its 
congruence subgroups in cases with discrete Wilson lines Ij^l ) , under which the matter fields transform according 
to their modular weights. For the Uk untwisted-sector states, these weights are ti = — 1 (if i = k) and (if i ^ k) for 
the three Kahler moduli. For twisted-sector states, ti = — f + Ti — Ni + Ni , where ri,Ni, Ni are the H-momenta and 
oscillator numbers. The superpotential has a weight ( — 1, —1, —1) (since the combination of Kahler potential and 
superpotential % + log |W[ 2 has weight (0,0,0)), so the modulus-dependent Yukawa couplings must have certain 
definite weights under the duality group and they are the modular forms. Physically, the modulus dependence 
accounts for the world-sheet instanton effects |k| and may be important for suppressing some Yukawa couplings. 
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for quantum numbers in the second and third complex planes. We note the above rules can be 
understood as discrete global symmetries (Z3 and Z2 respectively) in the field theory language. 

The conditions in the first plane are more complex since the fixed points need to be reorganized in 
terms of the #-eigenstates. Denoting the states by {^') a {i) , in addition to the apparent requirement 
that the product of ^-eigenvalues 



there are further conditions arising from eq. IC7I We list them as follows: 

1. There is no additional selection rules if the couplings involve states in the T\ sector - all 
gauge-invariant couplings consistent with eas. IC51 IC6I and IC8I are allowed. This follows from 
the fact that in the T\ sector all G2 root vectors are in the same conjugacy class. 

2. For couplings not involving any T\ twisted state, the selection rule is determined by a Z3 
(Z2 x Z2) discrete symmetry for the (T3) states. The T2 j 4-sector states are related to the 
Z 3 classes by (7 = l)i = [0], (7 = 1) 2 = [1] + [2] and (7 = -1) = [1] - [2], and the T 3 -sector 
states are related to the Z 2 x Z 2 classes by (7 = l)i = [00], (7 = 1) 2 = [10] + [01] + [11] and 
(7 = uj,uj 2 ) = [10]+7[01]+7 2 [11]. The allowed couplings can then be worked out by applying 
the appropriate multiplication tables of the conjugacy classes, [a] [b] = [(a + b) mod 3] and 
[0161] [0262] = [(01 + 02) mod 2, (61 + 62) mod 2]. 

The above selection rules can be used to determine the non-trivial Yukawa couplings. For 
example, consider three-point couplings (which correspond to renormalizable terms in the super- 
potential). From the H-momentum conservation eq. 1(^51 and the point group rule eq. 1(^61 we find the 
following set of allowed couplings in the 7L§ model (all of these states have zero oscillator numbers), 



Taking into account the space group rules, we find the complete list of three-point couplings in 
model Al involving fields with non-trivial representations under the SU4C factor of the PS group 
(we have used the field notation of table [HJ , 



(C9) 



[/1W3, T X T 2 T^ 21TXT4, T 3 T 3 U 2 , T 2 T A U 3 . 



(C10) 




^ (£4)0/3 X% , Yj ( C 4)a /A /B , Y ( C ' 4 )a /A f% i 





48 



(S»)a (qi)A (5l)B , Sl (C3)Aa(C 3 )^ , (Cll) 

a=0,l a+/3=0 mod 2 

where the "family" indices A, B = 1, 2 satisfy a Z2 condition, A + B = 0mod2. 
2. Allowed Yukawa couplings in the Al model 

In this subsection, we list allowed Yukawa couplings in the Al model, using the string selection 
rules of appendix IC II These operators all have the form 0S n , where are composite singlet 
operators that involve non-trivial PS fields, and S are the singlet fields. (Of course, one can also 
consider operators involving other observable- and hidden-sector fields.) We work out some of the 
lowest order allowed operators. We shall also follow the field naming scheme in table ITTl 

First notice certain structures exist for these Yukawa couplings. A composite operator of singlet 
fields is equivalent to the unit operator 1 if it is inert under all the global and local symmetries of 
appendix IC II Thus the superpotential factorizes W = {1}W, with the first few terms of {1} given 
by 

{1} = 1 + S19S32 + Sg(S4S32 + S5S33 + S10S26) + SfgS% 2 + ' - - 

= (1 - Si^r 1 !! - S 9 (S 4 S 32 + S 5 S 33 + SioSae)]- 1 x (• • •), (C12) 

where we have performed a "resummation" in the second line. Note that {1} can be absorbed into 
the Kahler potential by a transformation, DC — > 3<C + log |{1}| 2 . As a result, the sufficient conditions 
for F-flatness and vanishing cosmological constant are dW = W = 0. In the following expressions, 
we shall omit the caret. 

We list the relevant effective operators as follows, suppressing the A and a indices (defined in 
the last section) except for those of the matter fields and the x°i X° Higgs fields. 

A. Color triplet masses 

• The (6, 1, 1)(6, 1, 1) operators: 

w = (Si + s l9 s 25 + s 9 §f ] + S2S 2 a§ { 8 ) )(c 3 ) 2 + (S2S21 + s[ 2) sj, 2) )(c 4 ) 2 

+ (Sl^ + S^SgS^ + Si§>^ + Sig^S^) C3C4 

+ 5259522C2C4 + S^S^Sg 2 C2C3 + Sio§^\Ci) 2 + S10S4 2 C\C\ 

+ [Sio§> 3 2 ^ + Si 3 §^ + S1S10S12S32 + SgSio§^]CiC 3 + ■ ■ ■ (C13) 

• The (6, 1, 1)(4, 1, 2)(4, 1, 2) operators: 

W = C 1 ^ + (5 2 S 12 5 24 532 + 4 2) §4 2) )Ci(x?) 2 
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+ Sio(Sio§3 + <Si3§2 2 )Ci(X2) 2 + S^SgS^S^C^Xi) 2 

+ 5 26 (SS 2) + S 2 S 9 S 24 + Si§f)<7 3 (xD 2 + (Mo + S 10 S l9 S 25 + S 9 S l0 sf ] + SiS 2 3) )C3(X2) 2 
+ S 2 S^S 7& C^f + (5io§i 2) + S 6 S 7 S U S 18 + 5i5io4 2) + S 2 S? ) )C A (x c 2 ) 2 + ■ • • , (C14) 

• The (6, 1, 1)(4, 1, 2)(4, 1, 2) operators: 

w = c 3 r A x c a + c 4 (ftx c a + f c Af c B ) 

+ (S 9 + S' 1 5io5 21 5 22 + Sf ^[dfm + C 3 (/ 3 c Xa + /1/b)] 

+ (4 2) + s 10 (5 3 o4 2) + 5i 3 4 2) ))P 3 / 3 7l + C 4 (/ 3 C ) 2 ] + (4 2) + 5 26 §f)[C 3 x^^ + C A f A xll 
+ 5i 5i25 32 Ci/|xa + [5 9 S5 2) + 5i 3 (5i §i 2) + 5i 3 §( 2) )]C 3 (/ 3 c ) 2 

+ S w Si^S 2 iS 22 Cif§f A + S 9 Si Si 2 S 32 Ci(f§Xa + IaIb) + (S^^^XaX^ (C15) 
W( § ) d Si 2 (5 11< S 30 4 2) +5'io5 32 4 2) )C lX c a x c p + S 2 S 19 S 21 S 2 22 C 2X c a x C f3- (C16) 

• The (6, 1, 1)(4, 2, 1)(4, 2, 1) operators: 

W = S 2 SgSioSi 2 Ss 2 Cif3 + S 2 S 21 §>^C 2 f3 + S 2 SgCsf3 + S 2 C4f 2 

+ 5 2 5gSioS'l 2 S' 32 Ci/ 3 ,/U + S' 21 § 7 2 C 2 / 3 /^ + («S , 2 5g + S 2 lS23^)C3f3fA + S 2 C 4 fsfA 

+ SgSi(jSi 2 S3 2 CifAfB + 5'io(5'25'i2S , l 3 S'22'S'32 + 5' 3 5'^ 2 52 3 5 3 2 + »S'65'l4»S'l6'S'22'S' 33 )C2/ J 4/s 

+ (Sg + 5^S 2 2) + 5i5io52 1 52 2 )C 3 / j i/b + C 4 /a/b + • • • , (C17) 

In the above equations we have used the following composite singlet operators of dimensions 
two and three, 

Sj 2 = S4.S25 + S5S29, §2^ = ^45 32 + S5S33 + SiqS 2 q, 

Sg 2 "* = S12S25 + SisS 2 q + 5i75 32 + 5ig5 3 o, S 4 2 ^ = S4S30 + S5S31, 

S 5 2 = S4S1S + SiqSu, §^ = S 2 iS 2 /i + S 22 S 2 s, 

= S3S4 + S5S7, Sg 2 ' 1 = 5g + 5 2 l5 (C18) 
Sj 3 = S9S21S22 + S'io<S'i 2 5 3 o, § 2 3 ^ = >S'io»S'ig5' 3 2 + S'ii5i75 3 o, 
4 3 ' ) = 5'io5'2o»S' 33 + S'iiiS'i8<S' 3 i, sjj 3 ^ = 5 3 5i 3 5i7 + SgSxoS 2 4- (C19) 

-B. Quark and lepton Yukawa couplings 

Quark and lepton Yukawa couplings are operators of the (4, 2, 1)(4, 1, 2)(1, 2, 2) type. They 
are given by 

(a ,6) = fahifS, a, 6 = 1,2,3. (C20) 
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In addition, the PS symmetry breaking fields may enter the effective fermion mass operators at 
higher order. We thus consider the operators of the (4,1,2) (4,1,2) type given by 

Oi = x!x c a , 02 = %X c a (C21) 

Using these operators, we search for all the allowed higher dimensional operators of the form, 
{atb) ((Di) n S n ' , i.e., (4,2,1)(4, 1, 2)(1, 2, 2)[(4, 1, 2)(4, 2, l)] n S n ', with the smallest possible n, n' 
(for n + w'<4orn'<5). The result is (A, B = 1, 2) 

TV = §2 2 ' ) 0(3 ,A) + (<S'io5'l2S'24S'30 + ^9SwS22§>q'^ + 5*21 5*22 + S3S12O1 02)0(^4 3) 

+ (Sio^sf + S 3 S 9 S 12 Q2)Q(A,B)- (C22) 

C. Neutrino masses 

Consider the three PS breaking operators 

0y = X?3& i,J = l,2. (C23) 

We then find higher dimension effective Majorana neutrino mass operators of the form (with 
minimal number of singlets) 

W = (Sf 5 26 + 5 1 S 2 i52 2 )/3 C / 3 c 11 + S 13 / 3 7 3 c 0i2 + 5 9 5 10 8f /3 C / 3 C 22 + S 9 S 26 f-!f A 0n 

+ (5 10 4 2) + 5i 3 4 2) )/3 C /l0i2 + ^26/l/^0n + S 10 8i 2) /l/BOi2 + S 9 S w f c A f c B (D 2 2. (C24) 



APPENDIX D: GAUGE COUPLING EVOLUTION 



In this appendix, we derive the GQW equations for gauge coupling unification in the 7L§ string 
models. We shall work in the orbifold GUT limit and also make some simplification assumptions 
of the matter spectra. 

It is well known that gauge coupling unification in heterotic string models has a serious problem 
in the perturbative regime js^llssl]. From a simple reduction of the lOd string effective action, one 
finds G N = e 2 V 4 /128vrU, 

^string — e 2< ^ > Q; / ^/167rU, where <p is the dilaton, V the volume of the 
compactified space, and a' the Regge slope, thus 19 



Mpi = J Mstring. (Dl) 

y ^string 



A more careful definition, taking into account of the renormalization scheme dependence, is given in ref. |31 
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Assuming V ~ ^ G UT (where Mqtjt - 3x 10 16 GeV), a s trin g = «gut — V 24 an ^ ^ < 1, on e 
gets Mpj < Q!qut M Gut ^ 2 x 10 17 GeV, which is apparently incorrect. Note that this result relies 
on the assumptions that the compactified space is symmetric, i.e., all the dimensions are of similar 
sizes, and there are no additional states near the GUT scale. Both assumptions are invalid in our 
models. 

In heterotic models, it is well known that gauge couplings can obtain potentially large threshold 
corrections only if the models contain N = 2 sub-sectors 
exactly this kind of property, since both the Z2 and Z3 sub-orbifold twists leave exactly one complex 
plane unrotated. It is, however, a very complicate matter to compute gauge threshold corrections 
in string theory in the presence of discrete Wilson lines. (The simplest way to compute these 
corrections is using the target space duality symmetry [5^ . however in the presence of discrete 
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361 ]. Our 7L§ orbifold models have 



Wilson lines the duality groups are broken to their discrete subgroups 
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simpler Z3 models, threshold corrections are only known numerically 

Fortunately we can use a much simpler field theoretical method to compute gauge threshold 
corrections in the orbifold GUT limit, i.e., when one of the SO4 dimensions is much larger than the 
string length scale and other compactified dimensions. 21 In this limit, all the winding modes and 



32, 



36j one would 



KK modes in small dimensions have string scale mass, and according to refs. 
expect the dominant contributions to the threshold corrections come from the KK modes of the 
large dimension. From now on, we shall assume this is correct and neglect all the contributions 
from states with string scale mass. The gauge coupling at the string scale, a s tring> imposes a 
boundary condition for the renormalization group equations at the cut-off scale, which is taken to 
be the string scale, M str i ng . 



We follow the field theoretical analysis in ref. [6l| (see also 62j). It has been shown there the 



correction to a generic gauge coupling due to a tower of KK states with masses Mkk — m/R is 

°-'< A > = ■-'<«> 7 (D2) 

where the integration is over the Schwinger parameter t, (Aq and A are the IR and UV cut-offs, and 
r = 7r/4 is a numerical factor. #3 is the Jacobi theta function, 63(f) = ^m=-oo e 17 ™ 2 *, representing 
the summation over KK states. 



20 Threshold corrections in models with continuous Wilson lines |63[ have been computed in refs. |33| . Our models 
with discrete Wilson lines, however, are not the limiting cases of those with continuous lines. 

21 Admittedly the field theoretical calculation suffers from the usual UV divergence. The result is sensitive to the 
cutoff scale and needs be dealt with caution. However, we do not expect the RG evolution of the difference of 
gauge couplings to be affected much by our field theoretical treatment. 
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In our models, there are several modifications in the calculation. Firstly there are four sets 
of KK towers, with mass M KK = m/R (for P = P' = +), (m + l)/R (for P = P' = -) and 
(m + 1/2) /R (for P = +, P' = — and P = —, P' = +), where m > 0. The summations over KK 
states give respectively \ ^^(it/irR 2 ) — l) for the first two cases and ^92(it/irR 2 ) for the last two 
(where 02{t) — Sm=-oo e 17r ( m+1 / 2 ) 2 *), and we have separated out the zero modes in the P = P' = + 
case. Secondly, the PS symmetry in our models must be further broken down to that of the SM. 
This breaking, in principle, can be induced by brane or bulk states. However, in sect. UVBl we have 
shown this breaking is more likely due to non-renormalizable couplings of the states in the T2 t 4 
sectors. Since these states are identified with bulk states in the orbifold GUT limit (c.f. sect. lIII A"|) . 
we shall assume, in what follows, the breaking of PS to the SM is bulk breaking. We further assume 
the breaking scale, Mpg, is smaller than or equal to the compactification scale, M c , so that we can 
neglect mass corrections to the massive KK states. 

Tracing the renormalization group evolution from low energy scales, we are first in the realm 
of the MSSM, and the beta function coefficients are fr¥ SSM = When we pass the 

PS breaking scale, Mpg, the beta function coefficients become (&++ + &brane)i> where the two 
terms represent contributions from the bulk and brane states. These coefficients can contain 
contributions from additional states besides those of the MSSM, for example, from a vector-like 
pair (4, 1, 1) + (4, 1, 1). The next energy threshold is the compactification scale M c . From this 
scale to the string scale, we have the four sets of KK states. Since we have assumed Mpg <C M c , 
the beta function coefficients are those of the complete PS representations. 

Collecting these facts, and using # 2 (i*A-R 2 ) — ^zijt/^R 2 ) ~ sfjR for t/R 2 <C 1, we find the 
GQW equations, 

. . 27T i ^MSSM i MPS /,PS . i \ i Mitring 



+ 6 f ^ M log ^ + (&£ + 6 brane ), log ■ 



<Xi(fi) ^string A 4 M PS 

~ + ^ S _), log ^ + 6 S (^jgS - l) , (D3) 

for i = 1,2,3, where we have taken the cut-off scales, fio = M c and A = M s t r i ng . = 
^2p=± p'=± bp P ,, so in fact it is the beta function coefficient of the orbifold GUT gauge group, S- 
The beta function coefficients in the last two terms have an 3sf = 2 nature, since the massive KK 
states enjoy a large super symmetry. For S = SOio and Eg, & s = — 16 + 4ni6 + 2riio and — 24 + 6?t,27, 
where nio, n^Q and are numbers of bulk hypermultiplets in the respective representations. To 
accomplish bulk breaking, we need n-n = 4 in models Al/B and ni6 = 4, 2 < ri\Q < 6 in model 
A2, therefore 6 Ee = and 6 SO(10) = 2n 10 . Coincidentally there is no power-law running in the Eg 
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model. 
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